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FORTRAN PROGRAM TO GENERATE ENGINE INLET FLOW 
CONTOUR MAPS AND DISTORTION PARAMETERS 
by John H. Dicus 
Lewis Research Center 

SUMMARY 

A computer program has been written to generate gas turbine engine inlet flow 
contour maps and inlet flow distortion parameters based on any array of measurements 
describing the flow conditions at the engine inlet. 

As well as being able to choose the type of contour map that is generated, the user 
has the ability to construct any distortion parameter desired, as long as it depends 
solely on the conditions at the engine inlet. To accomplish this, a simple subroutine is 
written which combines the information generated by the program into a parameter. 

The report includes the FORTRAN IV computer program, a discussion of its 
operation and limitations, and two example cases with the input and output shown. 

A CFACE source deck on tape is available from COSMIC (Computer Software 
Management and Information Center), University of Georgia, 30601. The program can 
by ordered by using the number of this report as identification. 


INTRODUCTION 

This report presents a computer program to generate inlet contour maps and inlet 
distortion parameters for gas turbine engines. Inlet flow distortion poses a problem for 
all types of gas turbine installations: commercial, civil, and military. The compressor 
system of a gas turbine engine may stall when subjected to a distorted inlet flow field, 
whether it be a pressure distortion, a temperature distortion, or a distortion of some 
other flow parameter. Also, engine and airframe manufacturers need a method to 
insure compatibility between the inlet and engine of an aircraft before they are tested 
as a unit . 

Currently an intensive effort is being made to correlate a gas turbine engine's 
susceptibility to stall with a description of the flow field at the engine inlet. One approach 


to the problem is to generate inlet distortion parameters based solely on the measured 
flow conditions at the inlet to the engine and to attempt to correlate these distortion 
parameters with some measure of the loss in the engine stall margin. This procedure 
requires an accurate definition of the flow field at the engine inlet. Reference 1 presents 
the stall margin of an engine as a function of an inlet distortion parameter. In refer- 
ence 2, the value of a distortion parameter causing engine stall is shown as a function 
of corrected compressor massflow. 

The usefulness of the program presented in this report lies in the fact that it is 
sufficiently general to accommodate almost any inlet probe geometry. The only stipu- 
lation is that the probes must lie on some definable rings. The rings may vary in 
number and in radial location. They may contain different numbers of probes, and the 
probes may be located at any circumferential position on the rings . 

Included in this report is the FORTRAN IV program CFACE, complete with two 
example cases showing both input and output. Also included is a brief discussion of the 
methods used to generate the contour maps . Four examples are given to show how to 
build specific distortion parameters within a special subprogram named DISPAR. 

A CFACE source deck on tape is available from COSMIC (Computer Software 
Management and Information Center), University of Georgia, 30601. The program can 
be ordered by using the number of this report as identification. 


GENERAL PROGRAM DESCRIPTION 

The purpose of the program is to generate contour maps and distortion parameters 
that describe the flow at the measuring station immediately ahead of a gas turbine engine. 
However, it may be used to generate contour maps for any type of measuring station. 
Although total and static pressures and temperatures are the most common inputs to the 
program, provision has been made for the input of special measurements such as flow 
angles and turbulence levels. The FORTRAN IV program CFACE operates on these 
input measurements and the probe array geometry, which is also part of the input, to 
generate the contour maps and distortion parameters. The entire program listing and 
error messages are contained in appendix A. 

Also included in appendix A is a brief description of each subroutine telling what 
each one accomplishes. Also, comment statements are placed throughout the program. 
These two features enable the reader to follow the flow of the program and understand 
what is taking place in each section. 

CFACE generates inlet contour maps with a choice of mapping parameters. To 
provide contour maps and distortion parameters, interpolation between input measure- 
ment locations is done by a technique involving linear interpolation. The method of 
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interpolation is discussed in appendix B. The contour maps are represented by symbols 
(numbers and letters) on a picture produced by a line printer. Twenty mapping symbols 
are used and they are described by a key which accompanies each contour map. Interpre- 
tation of the contour maps is discussed in the section entitled OUTPUT . 

Three different types of contour maps are available when the input measurements 
are pressures. They are (P t , max, f “ P i> /P t, max, f ’ ^percent; (P. - P t , av , f)/ 

P. _ T . , in percent; and Mach number. (Symbol definitions appear in appendix C.) 

X j 3-V j X 

When temperature is the input, a map of (T. - T t av f )/T t ay ^ , in percent, can be 
generated. 

As mentioned before, provision has been made for the input of special measurements 
such as flow angles and (AP) /P. Calling this special variable "X", either 

i III o 

(X. - X f )/X in percent, or X ± n(AX) may be mapped. In the case of 
X ± n(AX), AX is specified by the user; X is the mean value of X rounded off to the 
nearest AX; and n = 0, 1, 2, ... , 10. If a map of pressure in P + n(AP) or 
temperature in T ±n(AT) is desired, the special input format can be used. 

CFACE generates a variety of simple circumferential and radial distortion 
parameters that enable the calculation of almost any specific distortion parameter as 
long as its calculation is based solely on the conditions at the engine inlet . To imple- 
ment this procedure, user-written FORTRAN IV instructions are inserted into a 
special subprogram named DISPAR. These instructions combine the general param- 
eters calculated by the program, the probe array geometry (part of the input), and any 
information unqiue to the specific distortion parameter (inserted by the user) into the 
desired form. Distortion parameters can be calculated regardless of the type of input. 
Although the user has access to the input data (in the common block DATTRN) within 
the subroutine DISPAR, care should be taken not to alter their values as they will be 
used elsewhere in the program. 

The output of all of the distortion parameters is handled by the output subroutine 
LNPOUT. A special section in LNPOUT is reserved for the output of the specific 
parameters. Another special section is reserved in LNPOUT for the corresponding 
output format statements. A special common block named PARTRN is reserved for 
transmitting the specific distortion parameters from DISPAR to LNPOUT. 

The simple distortion parameters calculated by the program are mentioned in the 
section entitled OUTPUT and are described in appendix D. Appendix D also contains 
four examples of user-defined specific distortion parameters and the FORTRAN IV 
instructions to calculate them. A listing showing the use of the common block PARTRN 
and the subroutines DISPAR and LNPOUT is shown in appendix A. 
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INPUT 


Two example cases showing input and output are contained in this report. They 
appear in appendix E . Appendix B shows the method of interpolation used to generate 
the contour maps. Referring to appendixes B and E while reading the following section 
should aid in its under standing. This particularly applies to figure 2, which is refer- 
enced in appendix B . 


Input Variables 


The input variables are shown below. A more detailed description of these variables 
is contained in the section entitled Instructions For Preparing Input. 

QQ (integer) indicates how many complete sets of measurements are to be 

read in 


WHICH (integer) 


RANGE (real) 
LPCENT (logical) 


DELX (real) 
STATI (locigal) 


controls the type of contour map to be generated 

WHICH = 1, map of (P tmaX;I - ^ t ,uax,l’ percent 
WHICH = 2, map of (Pj - P^ ay , f )/P^ ay { , percent 
WHICH = 3, map of Mach number 
WHICH = 4, map of (T. - T t f )/T t ay? f , percent 
WHICH = 5, map showing distribution of special variable - 
can be represented as (X i - X ay { )/X f (in percent) 
or as X + n(AX) 

total range of percentage represented on contour map when 
WHICH = 1, 2, or 4 or when WHICH = 5 and LPCENT = .TRUE.. 

controls whether (X A - X &y { )/X ay f or Xi n(AX) will be 
mapped when WHICH = 5 

LPCENT = .TRUE., causes (X^ - X fly f)/X &y f to be 
mapped 

LPCENT = .FALSE., causes X + n(AX) to be mapped 

magnitude of AX that each of the 20 mapping symbols will 
represent when WHICH =5 and LPCENT = .FALSE. 

indicates whether or not the innermost ring is treated as ring 
of static measurements 
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STATO (logical) 

LCNTFL (logical) 

RMACH (real) 

LDPAR (logical) 

NR (integer) 

RADLOC (10) (real) 
NP (10) (integer) 
ANGLOC (10, 20) (real) 

TITLE 1 (A -format) 
TITLE 2(A -format) 

P (10, 20) (real) 


STATE = . TRUE . innermost ring is ring of static measurements 
STATI = .FALSE, innermost ring is not ring of static measure- 
ments 

indicates whether or not the outermost ring is treated as a ring 
of static measurements 

STATO = .TRUE, outermost ring is ring of static measure- 
ments 

STATO =-. FALSE, outermost- ring is not ring of static measure- 
ments 

controls whether or not center of contour map will be filled in 
LCNTFL = .TRUE, center of map will be filled in 
LCNTFL = .FALSE, center of map will not be filled in 
(interpolation will stop with the innermost ring) 

maximum Mach number that will be represented on contour map 
when WHICH = 3 

controls whether or not distortion parameters are calculated 
from the input measurements 

LDPAR = .TRUE, distortion parameters are to be generated 
LDPAR = . FALSE . distortion parameters are not to be 
generated 

number of rings of measurements to be read in (10 is the maxi- 
mum), numbered from innermost ring out 

array of radii, one for each of the NR rings, cm (in. ) 

array containing number of measurements on each of the NR rings 

array containing angular position (in deg) of each measurement 
, on each of the NR rings (There may be up to 20 angular posi- 
tions on each ring.) 

can contain up to 80 alphanumeric characters for general title 
information 

can contain up to 80 additional alphanumeric characters for 
general title information 

array containing value of each input measurement on each of NR 
rings (There may be up to 20 measurements on each ring.) 
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Instructions For Preparing Input 


QQ, WHICH, RANGE, LPCENT, DELX, STATI, STATO, LCNTFL, RMACH, and 
LDPAR are read in under the namelist INIT and go on the first two cards. These two 
cards should not be repeated for each set of input measurements; they only appear once. 

QQ. - QQ may assume any positive integer value and tells the program how many 
contour maps are to be generated. 

WHICH, RANGE, LPCENT, DELX - A value for RANGE is needed when WHICH = 1, 
2, or 4 or when WHICH = 5 and LPCENT = .TRUE. . When WHICH - 3 or when 
WHICH = 5 and LPCENT = .FALSE., RANGE may be omitted from the namelist. There 
are 20 mapping symbols, and the value of RANGE divided by 20 will be the spread in 
percentage represented by each mapping symbol. RANGE =20.0 will cause each map- 
ping symbol to represent a spread of 1 percent. (When WHICH = 1, a total range of 
20 percent would correspond to 0 to -20 percent. In all other instances, a total range 
of 20 percent corresponds to -10 to +10 percent.) This is usually sufficient for most 
pressure or temperature distributions. If RANGE equals 30.0, each mapping symbol 
will represent a spread of 1.5 percent. Obviously, care should be taken to avoid 
choosing values of RANGE that would cause each mapping symbol to represent an in- 
convenient spread of percentage . 

When WHICH = 5, it is necessary to set LPCENT equal to either .TRUE, or 
. FALSE . . If WHICH * 5, LPCENT can be omitted from the namelist as its value is 
unimportant. If WHICH = 5 and LPCENT = . FALSE . , DELX must also be given a 
value (otherwise DELX can be omitted from the namelist) . Since there are 20 mapping 
symbols, DELX should be chosen such that 20 times DELX is greater than or equal to 
the range of measurements being mapped. Consideration should also be given to pick- 
ing a reasonable DELX such as 0.02 or 30.0, as opposed to 0.0218 or 31.66. 

STATI, STATO - When the innermost ring is a ring of static measurements, STATI 
should be set equal to .TRUE. . Then the measurements on the innermost ring will be 
excluded in the calculation of the average and in the search for the minimum and the 
maximum. They will also be treated as static measurements in the calculation of 
distortion parameters . In general, a ring of static measurements will be treated as 
just another ring of measurements in all other respects (such as interpolation) . The 
same comments apply for the variable STATO if the outermost ring is a ring of static 
measurements. Since the calculation of Mach number requires static pressures, 

WHICH = 3 cannot be chosen when both STATI and STATO are set equal to . FALSE . . 

LCNTFL - When the innermost ring is not a ring of static values, the user may 
want the mapping to continue through the center of the profile. If this is the case, 
LCNTFL should be set equal to .TRUE. . 
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If the user wants the mapping to stop with the innermost ring, LCNTFL should be 
set equal to .FALSE. . The condition of LCNTFL and STATI both equal to .TRUE, is 
a condition that cannot be allowed. 

RMACH - RMACH determines the maximum Mach number that will be represented 
on the Mach number contour map when WHICH = 3 . Each of the 20 mapping symbols 
will represent a delta Mach number equal to RMACH -f 20. Care should be taken to 
select RMACH large enough, but of such a value that the delta Mach number will not be 
awkward. If WHICH * 3, RMACH can be omitted from the namelist. 

LDPAR - If the calculation of one or more specific distortion parameters is 
desired or even if only the calculation of the simple distortion parameters is desired, 
LDPAR should be set equal to .TRUE. . Otherwise, LDPAR should be set equal to 
.FALSE. . The calculation of distortion parameters is usally desired when the input 
measurements are pressures. However, the situation may occur where distortion 
parameters are desired when the input measurements are temperatures or some other 
measurement. 

The input of NR, RADLOC, NP, and ANGLOC will be subject to specific input 
formats and should all be on separate cards. They should not be repeated for each set 
of measurements; they should appear only once. 

NR - The number of rings of measurements should be read in as NR with an (12) 
format. Rings of static measurements should be counted. The maximum number of 
rings allowed is 10. 

RADLOC - The radius of each ring goes into the array RADLOC. The radii may be 
in either centimeters or inches. RADLOC(l) corresponds to the innermost ring while 
RADLOC(NR) corresponds to the outermost ring. The input format is (F5. 2,9 (1X,F5. 2)) . 

NP - The array NP contains the number of measurements on each ring. Again, 
NP(1) contains the number of measurements on the innermost ring. The maximum 
number of measurements that each ring can contain is 20. There may be as few as 
one measurement on a ring. The input format for NP is (12, 9(1X,I2)). 

ANGLOC - The angular position of each measurement on each ring goes into the 
array ANGLOC. The first subscript pertains to the ring number, and the second 
subscript pertains to the measurement referred to on that ring. For example, 
ANGLOC(5, 14) stores the angular position of the 14th measurement on the 5th ring. 

The values in ANGLOC should be in degrees, and there can be up to 20 values on each 
ring. The input format for ANGLOC is (F5. 1, 9(1X, F5. 1)), hence, only 10 values can 
go on a card. Each ring has its own card (or cards if there are more than 10 values 
on a ring) associated with it. Angular position is measured clockwise from the 
12 o'clock position. On a given ring, the first measurement falling on or after the 
12 o'clock position is the first measurement on that ring (see figure 2 in appendix B). 

TITLE 1, TITLE2, and P also go on separate cards. However, they appear more 
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than once. The number of sets of these variables should be equal to the number of maps 
to be generated (QQ) . 

TITLE 1, TITLE2 - TITLE 1 is a string of alphanumeric characters which is read in 
under an (80A1) format. The title can be used for general information such as airflow 
or reading number. This information will appear as part of the output. The purpose of 
TITLE2 is the same as TITLE1. It also has an (80A1) format. 

P^ - The measurements that are to be mapped should be read into the array P. As 
in the case of the subscripts of ANGLOC(5, 14), P(5, 14) contains the 14th measurement 
on the 5th ring. There can be up to 20 measurements on each ring. The contents of 
ANGLOC correspond to the contents of P. For example, ANGLOC(5, 14) is the angular 
position of P(5, 14). 

The input format for P depends on the value of WHICH. When WHICH = 1, 2, or 3, 
the input format for P is (F7. 3, 9(1X, F7. 3)), where P is expressed in kilopascals (kPa) 
or psi. When WHICH = 4,. the format is (F7.2, 9(1X, F7.2)), where P is expressed in 
K or °R. And when WHICH = 5, the format is (F11.5,4(1X, F11.5), where the choice of 
dimensions for P is left to the user . Except for when WHICH = 5, the values of up to 
10 measurements can be placed on one card. If there are more than 10 measurements 
on a ring, a second card for that ring will be used. Therefore, each ring will have one 
or two cards corresponding to it. If WHICH = 5, the values of only five measurements 
will fit on one card, and up to four cards for a given ring may be used. 


OUTPUT 

Two examples of program output are shown in appendix E . The first page of the 
output contains input information, specifically, the value of each measurement, its 
angular position, and the value of the mapping parameter corresponding to each 
measurement. The first section of this page lists the angular position of each measure- 
ment for each ring radius. Section two lists each measurement on each ring. There is 
a one-to-one correspondence between the angular positions in the first section and the 
measurements in the second section. The third section shows the value of the parameter 
being mapped for each measurement location. Again, there is a one-to-one correspond- 
ence between the values in this section and the values in the other two sections. 

Page two of the output contains the contour map as well as a key to the mapping 
symbols. A line can be drawn (by hand) between two adjacent regions having different 
mapping symbols. This line is a line of constant mapping parameter which can be 
determined from the key. The contour map is to scale with 19 line printer rows (or 
lines) equalling the radius of the. outermost ring of input measurements. 

The next two pages (pp. 3 and 4) of the output exist only if distortion parameters are 
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desired (only if the logical variable LDPAR has been set equal to .TRUE.). The first 
section on page three of the output lists eight simple parameters dealing with the cir- 
cumferential distortion. They are listed for each of the 10 rings. Section two of this 
page lists six simple parameters dealing with radial distortion. They are listed for 
eight discrete circumferential locations. The method used to enable calculation of these 
simple radial distortion parameters is discussed in appendix F. On this same page there 
are also seven simple distortion parameters pertaining to the complete array of input 
measurements. Appendix D discusses the calculation of the simple distortion param- 
eters. — 

The last page of the output (p. 4) contains information dealing with the specific 
distortion parameters that can be defined by the user. The format of this page is also 
user-defined. An example of the output is shown in appendix E, and a sample listing is 
shown in appendix A. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 3, 1973, 

501-24. 
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APPENDIX A 


PROGRAM LISTING AND ERROR MESSAGES 

A block diagram showing the program flow through the various subroutines is shown 
in figure 1 . 



Figure 1. - Calling relation of subroutines. 


Once the program flow reaches CALL PF1X(. FALSE.) in CFACE, and if 
LCNTFL=.TRUE. , all input data are moved outward by one ring, and a dummy ring 
takes the place of ring one . This happens so that the program can interpolate through 
the center of the contour map. For example, when the program flow reaches CALL 
PFIX(. FALSE.) in CFACE (if LCNTFL=.TRUE.), ring NR becomes ring NR+1, ring 1 
becomes ring 2, and so on. Also in this call to PFIX, the input measurements are 
transformed into mapping parameters . Therefore, the user must take care not to 
access the input data in any subroutine other than DISPAR. 


CFACE 


CFACE is the main driving program. It controls the number of cycles through the 
program and calls the subroutines in the proper order. 
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MAIN DRIVING PROGRAM 
COMMON/CYCLE/LREAD 

COMMON/READTR/STATl, STATO, LCNTFL, LDPAR 

C OMMON/ WR ANG E/ X BAR, RM AV, WHICH, RANGE, LP CENT, RMACH, HE LX 

COMMON/ F I N I /LQU I T 

INTEGER WHICH 

INTEGER QQ 

LOGICAL STATI .STATO, LCNTFL, LDPAR 
LOGICAL LREAD 

LOGICAL LPCENT-.LOUIT ------ - _ 

NAMELI ST/ INI T/QQ, WHICH, RANGE, LPCENT, DEL X, ST AT I ,ST ATO.LCNTFL, RMACH 
A LDPAR 
LRE AD= .TRUE . 

LQUIT=. FALSE. 

READ IN CONTROL PARAMETERS 

READ ( 5 , I N I T ) 

DO 13 1=1,00 

CALL SUBROUTINE TO READ IN GEOMETRY AND MEASUREMENTS 
FOR THE INPUT MEASUREMENT PLANE 

CALL CARDRD 
IF(LQUIT) GO TO 10 

CALL SUBROUTINES TO CALCULATE DISTORTION PARAMETERS 

IF(LDPAR) CALL PFIXI.TRUE.) 

I F ( LDPAR) CALL DISTRT 

CALL SUBROUTINE TO TRANSFORM INPUT MEASUREMENTS INTO 
MAPPING PARAMETERS 

CALL PFIXI .FALSE. ) 

I F ( LQU IT) GO TO 10 

CALL SUBROUTINE THAT CONTROLS INTERPOLATION PROCESS 
IN INPUT MEASUREMENT PLANE 

CALL MA I NLP 

CALL SUBROUTINE THAT CONTROLS THE OUTPUT 


CALL LNPOUT 
10 CONTINUE 
STOP 
END 


ooo noo onoooonoooo 


Subroutine CARDRD 


The reading of input data from cards is controlled by this subroutine. The output 
of probe angular locations and input measurement values is handled by CARDRD. 
"PROGRAM WILL NOT WORK. NR LESS THAN 2" 

will be printed out and the program will stop if there are less than two rings of input 
measurements . 


SUBROUTINE CARDRD 

SUBROUTINE TO READ IN GEOMETRY AND MEASUREMENTS FOR 
THE INPUT MEASUREMENT PLANE 

SUBROUTINE ALSO CONTROLS SOME OUTPUT 

COMMON/CYCLE/LREAD 

NP, RADLOC, ANGLOC, AND P HAVE THE FIRST 

ARRAY DIMENSION EQUAL TO 11. 10 OF THE 

11 POSITIONS ARE AVAILABLE FOR USE. THE 

REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 

COMMON/ DAT TRN/P! 11 ,20 ) , NR,R ADLOC! 111, ANGLOC ( 11,20), 

A NP<11) ,FACE< 39, 65> 

COMMON/RFADTR/ST AT I , ST ATO,LCNTFL,LDPAR 

COMMON/ WR ANGE/ X BAR, RM A V, WHICH, RANGE, L PC ENT, RMACH, DEL X 
COMMON/SCALF/FACTOR V 

COMMON/HEADG/TITLEl! 80),TITLF2( 80) , 

COMMON/FI NI/LQUIT 
INTEGER WHICH 
LOGICAL LREAD 

LOGICAL ST AT I ,STATO, LCNTFL ,LDPAR 
LOGICAL L PC ENT , LQU I T 

SECTION THAT READS IN THE GEOMETRY OF THE MEASUREMENT PLANE 

I F ( .NOT. LREAD) GO TO 20 
LREAD=. FALSE. 

RE AD ( 5 , 100 ) NR 
IFINR.GE.2) GO TO 31 
WRI TE ( 6,130 ) 

LQU I T= . TRUE • 

RETURN 
31 CONTINJE 

REAU(5,1D1) (RADLOCI I ) , 1=1, NR) 

F AC TOR=R ADLOC! NR 1/19.1 
RE AD (5, 102) (NP( I ) , 1=1, NR) 

DO 1C J=1 , NR 
K=NP( J) 

READ(5,103) ( ANGLOC! J, I ), 1 = 1, X ) 

10 CONTINJE 
20 CONTINJE 

SECTION THAT READS TITLES AND INPUT MEASUREMENTS 

* 

READ(5,12A) TITLE1 
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READ(5,124) r I TLE2 
DO 3D J=l,Nk 
K = NP ( J ) 

IF(WHi:H.NE.4.AND.WHICH.ME.5) READ(5,104) (P( J, I) , 1=1, K) 

I F ( WHICH. EQ. 5) READ(5,107) ( P ( J , I ) , [ = 1 , K ) 

IF(WHi:H.EQ.4) READ(5,106> (P(J,I)*I=1*K) 

30 CONTINUE 

WRI TE( 5,126 ) 

WRITE(S»125) TITLE1 
’ WR I T E ( 5 * 1 2 5 ) TITLE2 
WR I TE ( 5 , 1 09 ) 

SECTION THAT WRITES THE ANGULAR POSITIONS OF THE 
INPUT MEASUREMENTS 

WRITE! 5,121 ) 

DO 41 J = 1 , N R 
K = NP('J) 

IF(WHI:h.NE. 5) WRITE ( 6,1 23 ) RADLOC ( J ) t ( ANGLOC ( J , I) -, I =1 , K ) 

IF! WHICH. EQ. 5 ) WRITE(6,122) RADLOC < J >,( ANGLOC ( J , I ), I =1 , K ) 

41 CONTINJE 

WRI TE( 5,109) 

SECTION THAT WRITES THE INPUT MEASUREMENTS 

IF! WHI :H.NE.4.AND.WHICH.NE.5 ) WRITE! 6, 113) 

IF{WHICH.EQ.4) WRITE (6*114) 

I F ( WHI CH. EQ. 5 ) WR I TE ( 6 * 1 15 ) 

DO 40 J= 1 , NR 
K=NP(J) 

IF(WHI;h.NE.4.AND.WHICH.NE.5) WRITE! 6, 110) RADLOC! J),(P(J, I ) ,1=1, K 
A) 

IF(WHi:H.E0.4) WR I TE ( 6 , 1 1 1 ) RADLOC! J ) , ( P ( J, I ) , I =1 ,K) 

IF(WHi:H.EQ.5) WRITE (6,112) RADLOC! J),(P(J, I ), 1=1, K) 

40 CONTINUE 

WRI TE! 6,109 ) 

SECTION TO WRITE TITLES FOR NEXT SECTION OF OUTPUT 
NEXT SECTION OF OUTPUT IS CONTROLLED BV THE 
SUBROUTINE • PF I X ' 

IF(WHICH.EQ.l) WR I TE ( 6, 1 16 ) 

IF(WHICH.E0.2) WR I TE ( 6 , 1 17 ) 

I F ( WHI CH, EQ. 3 ) WRITE (6, 118) 

IF! WHICH. EQ. 4) WRI TE ( 6, 1 19 ) 

I F ( WHICH. EQ.5. AND. LPCE NT ) WRITE(6,120) 

100 FORMAT ( 12 ) 

101 FORMAT (F5.2 ,9(1X»F5.2) ) 

102 FORMAT! 12,9! IX, 12) ) 

103 F0RMAT(2(F5.1,9(1X,F5.1)/) ) 

104 FORMAT ( 2 ( F7 . 3 , 9 ( 1 X , F 7 . 3 ) / ) ) 

106 F0RMAT(2(F7.2,9( 1X,F7.2)/) ) 

107 FORMAT (4! FI 1.5, 4! IX, F 11. 5)/) ) 

109 FORMAT ( 1H0 ) 

110 FORMAT (5X,F5.2,4X, ID! 3X, F7.3)/14X, ID! 3X,F7.3) ) 

111 FORMAT (5X,F5.2»4X, ID! 3X , F7.2 1/14X, 10! 3X,F7.2) ) 

112 FORMAT (5X,F5.2, 4X , 5 ( 2X,F11.5l/3(l4X,5(2X,F11.5)/)> 

113 FORMAT ( IX/ 4X»oHRADIUS , 15X, BHPRESSURE/ ) 

114 FORMAT! IX/4X,6HRAD I JS, 15X, 1 1 HTEMPER ATURE / » 
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115 FORMAT (1X/AX,6HRADIJS* 15 X , 13HSPECIAL INPUT/) 

116 F OR MAT ( 1X/AX,6HRADIJS» 1 5X , 23HIPI-PMAX ) / PM AX » PERCENT/) 

117 F OR MAT ( 1X/AX,6HRADIJS» 1 5 X, 2 3H( PI-PAVG)/PAVG* PERCENT/) 

118 FORMAT ( 1 X/ AX , 6 HR AD I'JS»15X.»1 1HMACH NUMBER/) 

119 FORMAT ( 1X/AX,6HRADIJS» 15X» 23H( T I-TAVG)/TAVG, PERCENT/) 

120 F0RMAT(1X/4X,6HRA0IJS,15X,32HSPECIAL, ( X I -XAVG ) /XAVG » PERCENT/) 

121 FORMAT ( 1X//AX,6HRADI USi 1 5X , A3HANGUL AR POSITION OF INPUT MEASUREMEN 
ATS, DEG ) 

122 FORMAT (5X,F5.2»4X,5( 8X , F 5. 1 )/3( 14X»5(8X»F5.1)/ ) ) 

123 FORMAT (5X,F5.2*4X, 10 (5X, F5.1 ) / 14X , 1 0 ( 5X,F5.1 ) ) 

1 2 A FORMAT ( 80 A1 ) 

125 FORMAT (6X,80A1 ) 

126 FORMAT ( 1 H 1 ) 

130 FORMAT ( IX , 37HPROGRAM WILL NOT WORK. NR LESS THAN 2) 

RETURN 

END 


Subroutine DISTRT 


DISTRT generates all of the simple circumferential and radial distortion param- 
eters that are used in DISPAR. DISTRT calls INTERP and DISPAR, in that order. 


’’RING HAS ONLY ONE PROBE” 


will be printed as a reminder when the subroutine tries to calculate 0" for a ring which 
contains only one probe. 

”NO STATIC PRESSURES, NOT POSSIBLE TO CALCULATE” 

’’SPECIFIC DISTORTION PARAMETERS REQUIRING STATIC PRESSURE” 

will be written when DISTRT finds no rings of static pressure measurements as part of 
the input. The average static pressure SMAVG will be set equal to 5. 0 and all calcu- 
lations using SMAVG will be meaningless. 


SUBROUTINE DISTRT 
C 

C SUBROUTINE TO CALCULATE SIMPLE DISTORTION PARAMETERS 
C 

C - 

C NP, RADLOC, ANGLOC, AND P HAVE THE FIRST 
C ARRAY DIMENSION EQUAL TO 11. 10 OF THE 

C 11 POSITIONS ARE AVAILABLE FOR USE. THE 
C REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 

C 

COMMON/DATTRN/PI 11,20) , NR, R ADL OC ( 1 1 ) , ANGL DC ( 1 1 , 2 0 ) , 

A NP(ll) ,FACE(39,65) 

COMMUN/DI STRN/CMMINI ID ) , CMMAXf ID ) , CM A VG ( 1 D ) , CDPOP II 10 ) , CDP0P? 1 1 3 ) , 
A CDPOP3I 10) ,THTMIN( 10 ), THTPLSI 10) ,RMMIN( 8 ) , RMMA X ( 8 ) , RMAVGt 8 ) , 

B RDPOP1I8) , RDP0P2 ( 8) , RDP0P3I 8) , TMM I N , TMMAX, TMAVG, SMAVG, TDP1 , 
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C TDP2,T0P3,QL0C(10,20 > 

COMMON/RE ADTR/ ST ATI , STATCJ, LCNT FL , L DP AR 
COMMON/ l NTRN/R, THETA, PVALJE 
DIMENSION SEG( ID ) ,RADP< 8,10) 

DIMENSION THTTOT ( ID) 

LOGICAL LADD 

LOGICAL STATI ,STAT0, LCNTFL,LDPAR 
CALCULATION OF SIMPLE CIRCUMFERENTIAL DISTORTION PARAMETERS 


Nl = 1 
N2 = NR 

-I F ( STATI 1 Ni = 2 
IFISTATO) N2=NR-1 
DO 13 I =N 1 , N2 

CALCULATION OF CMMAX, CMMIN, AND CMAVG FOR EACH RING 

C MM AX ( I ) = 0. 0 
CMMINI I) =50000.0 
CMAVGI I 1=0. D 
JL=NP( I) 

DO 11 J=1,J1 

IF( P( I , J) . G T. CMMAX ( I ) ) CMMAX ( I ) =P( I , J ) 

I F ( P ( I , J ) .LT. CMMINI I ) ) CMMINI I )=P< I, J) 

CMAVGI I ) = CMAVG I I ) + ( P < I , J ) / FLOAT ! J 1 ) ) 

11 CONTINJE 

CALCULATION OF CDPOP'S FOP. EACH RING 

CDPOP1I I ) = ( CMMAX I I ) -CMMINI I ) I /CMMAX! I ) 

CDPOP2 I I ) = 1 CMMAX! I )- CMMINI I ) ) /CMAVG! I ) 

CDPOP3 I I ) = ( CMAVG! I ) -CMMIN! I ) ) /CMAVG! I ) 

10 CONTINUE 

CALCULATION OF TMMAX AND TMMIN 


TMMAX=3.3 
TMM I N= 50O0Q • 0 
DO 7 I =N1 , N2 

I F I CMMAX I I ) . GT. TMMAX ) TMMAX = CMMAX! I ) 
IFICMMINI I ) .LT.TMMIN) TMMI N = CMM I N ( I ) 
7 CONTINJE 

CALCULATION OF TMAVG 

D I V = 0 . 0 
PTSUM=0.0 
DO 33 I =N 1 , N2 
J 1 = N P I I ) 

DO 33 J=1,J1 
PTSUM=PTSUM^-PI I , J ) 

33 DI V = DI V+-1 . 0 

TMAVG=PTSUM/DIV 

CALCULATION OF TOP'S 

TDP1=( TMMAX-TMMIN) /TMMAX 
TDP2 = I TMMAX-TMMIN) /TMAVG 
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T0P3=( TMAVG-TMMIN J/TMAVG 

SECTION TO CALCULATE THTMIN FOR EACH RING 

DO 9 I = N 1 » N 2 
DO 13 LL= 1 f 1 0 

13 S EG I LL ) =0 . 0 
J 1 = NP ( I ) 

I F ( Jl.EQ.l) GO TO 8 
L = 1 

LADU=. FALSE. 

DO 12 J=1,J1 
KO= J-l 

IFIK0.EQ.0) K0 = J 1 
Kl = J 

A1 = ANGLOC (I * K 1 > 

K2= J + l 

IFIK2.GT.J1) K2=K2-J 1 
A2= ANGLOC I I ,K2) 

I F ( ( J*1).GT.J1) A2=A2+360.0 
IF(J.GT.l) GO TO 14 
I F ( P( I ,K1 ). LT.CMAVGI I ) ) GO TO 21 

IFIPI I, Ki).6Q. CMAVGI I ) • AND.P l 1 * K2 ) .LT.CMAVGI I ) ) GO TO 22 
GO TO 14 

22 CONTINUE 

I F ( P ( I , K0 ) . GE. CMAVGI I ) ) GO TO 14 
21 CONTINUE 
LADD=. TRUE. 

14 CONTINUE 

LOGIC TO DETERMINE EACH SEGMENT'S CONTRIBUTION TO THTMIN 
ON EACH RING 

R I NC=0 • 0 

I F ( P ( I . K2 ) .LT.CMAVGI I ) .AND. PI I , K1 ). LT.CMAVGI I ) ) RINC = A2-A1 
I F( PI I ,K2 ). LT.CMAVGI I ) . ANO . P I I , K.l ) . EQ. CMA VG I I ) ) R INC = A2-A1 
IF! PI I ,K2 J.EQ.CMAVGI I ) .AND. PI I,K1 ) .LT.CMAVGI I ) ) R INC = A2-A1 
IF! PI I ,K2).GT. CMAVGI I ) .AND.P I I, Kl ) .LT.CMAVGI I ) ) 

A RINC = ( I A2-AD/I PI I,K2)-P( I,K1) ))*l CMAVGI I )-P II ,K1 ) ) 

I FI PI I ,K2 ). LT.CMAVGI I ) .AND. PI I ,K1 ) .GT. CMAVGI I ) ) 

A R I NC= I (A2-AD/IPI I,K1)-PI I,K2) ) ) *1 CMAVGI I ) -P 1 1 , K2 ) ) 

LOGIC TO DETERMINE HOH TO ADD TOGETHER THE 
CONTRIBUTIONS OF EACH SEGMENT 

IF(RINC.EQ.O.O) GO TO 20 

IFIPII, K2J.LT. CMAVGI I). AND. PI I ,K1 ) .GT.CMAVGI I) ) GO TO 23 
24 CONTINUE 

SEGI L) =SEG( L)+RINC 
GO TO 25 

23 CONTINUE 

IFIPII ,K0).GE. CMAVGI I ) ) GO TO 24 
IFIJ.E9.1) GO TO 24 
L = L + 1 
GO TO 24 
20 CONTINUE, 

I F I J.EQ.l ) GO TO 25 

IF I SEGIL) .EQ.0.0) GO TO 25 

L =L+ I 
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25 CONTINUE 
12 CONTINJE 
K S - 

THTTOT ( I )=D.O 

00 18 K=1 » 1 0 

I F ( SEG ( K) .NE.0.0) KS=KS+1 
THTTOT !I> =THTT0T( I)*SEG!K) 

18 CONTINUE 

I F I LADD) SEG(1)=SEG( 1 )+SEG! KS) 

I F ( LADD) KS=KS-1 

SEARCH FOR LARGEST DEPRESSION BELOW 
AVG ON EACH RING 

THIS WILL BECOME THTMIN FOR 
THAT RING 

THT MINI I ) =0 . 0 
DO 19 K=1 , KS 

I F ! SEG < K) • GT .THTMINI I)) THTMINI I ) = SEG < K ) 

19 CONTINUE 
GO TO 9 

8 CONTINJE 
WRI T E ( 5 » 1 00 ) 

100 FORMAT ( 1X»25HA RING HAS ONLY ONE PROBE) 

THTMINI I >=3.0 

9 CONTINUE _ J 

CALCULATION OF THTPLS FOR EACH RING 

DO 6 I =N1 , N2 

THT P LSI I )=363.0-THTT0T( I ) 

6 CONTINJE 

CALCULATION OF AVERAGE STATIC PRESSURE 

DIV=0. D 
S TA TSM = 0. 0 

IF! . NOT. STAT I ) GO TO 29 
I 1=MPC 1 1 
DO 30 1 = 1,1 1 

30 STATSM=STATSM+P( 1, I) 

DIV=DIV+FLOAT I II ) 

29 CONTINJE 

I F ( . NOT .STATO ) GO TO 31 
1 1=NP< NR) 

DO 32 1=1,11 

32 STATSM = STATSM+P( NR, I ) 

DIV=DI V+FLOAT! II) 

31 CONTINJE 

IFIUIV.EQ.O. 0) DI V=1 . 0 
SMAVG=STATSM/DI V 
IF(SMAVG.EQ.O.O) SMAVG=5.0 
IF(SMAVG.EQ.D.O) WRITEI6,101) 

IFISMAVG.EO.O.O) WRITEI6,102) 

101 FORMAT! IX, AfaHNO STATIC PRESSURES, NOT POSSIBLE TO CALCULATE) 

102 FORMAT! IX, 56HSPECIFIC DISTORTION PARAMETERS REQUIRING STATIC PRESS 
AURE ) 

CALCULATION Of SIMPLE RADIAL DISTORTION PARAMETERS 
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DO 40 I R = 1 » NR 
DO 40 I THT= 1 , 8 
R=R ADLOCI IR ) 

TH£TA = ( I THT-1 ) <=45 

CALL INTERPOLATION SUBROUTINE TO OBTAIN INTERPOLATED VALUES 
OF MEASUREMENTS ON EACH RING AT DISCRETE THETA LOCATIONS 
NOTICE THAT THE SUBROUTINE * INTERP' HILL GO THROUGH 
THE PROCESS OF RADIAL INTERPOLATION EVEN THOUGH IT 
IS NOT NECESSARY. THIS REDUNDANCY IS NOT EXTREME 
SINCE THE RADIAL INTERPOLATION IS NOT EXECUTED 
TOO MANY TIMES. 

CALL INTERP 

RADP( I THT , IR ) = PV ALUE 

40 CONTINUE 

CALCULATION OF PMMAX , PMMIN , AND PMAVG FOR 
EACH DISCRETE THETA LOCATION 

DO 41 1=1,8 
RMMAXI I ) =0. 0 
RMM INI I 1=5000000.0 
RMAVGI I ) =0. 0 
DO 42 J= 1 , NR 

I F ( R ADP ( I , J ) .GT. RMMAXI I ) ) RMMAXI I )=RADP( I, J) 
IF(RADP(I,J).LT. RMM I Nil)) RMMINI I )=RADP( I , J ) 

RMAVGI I )=RMAVG( II+RADPt I,J l/FLOATI MR) 

42 CONTINUE 

CALCULATION OF RDPOP'S FOR EACH DISCRETE THETA LOCATION 

RDPOP1 I I ) = ( RMMAXI I ) -RMMINI I ) ) /RMMAXI I ) 

RDPOP2I I ) = ( RMMAXI I ) -RMM INI I ) ) /RMAVGI I > 

RDPOP3I I ) = ( RMAVGI I ) -RMM INI I ) ) /RMAVGI I ) 

41 CONTINUE 

CALL SUBROUTINE TO CALCULATE SPECIFIC DISTORTION PARAMETERS 

CALL DISPAR 

RETURN 

END 


Subroutine DISPAR 

This subroutine is the one which contains user-written program steps used to cal- 
culate specific distortion parameters. 

"IDC-MAX AND IDR-MAX CANNOT BE CALCULATED” 

"THERE ARE NOT EXACTLY 5 RINGS OF TOTAL PRESSURE MEASUREMENTS” 
will be written when the program attempts to calculate IDC^,,,^ and IDR^,^ when there 

niaX ni3x 
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are less than five rings of total pressure input measurements. 


SUBROUTINE DISPAR 
C 

C SUBROUTINE TO CALCULATE SPECIFIC DISTORTION PARAMETERS 
C 

c 

C NP, RADLOC, ANGLOC, AND P HAVE THE FIRST 
C ARRAY DIMENSION EQUAL TO 11. 10 OF THE 

C 11 POSITIONS ARE AVAILABLE FOR, USE. THE 
C REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 

C 

COMMON/DATTRN/PI 11,20) , NR, RADLOC ( 11), ANGLOC I 11,20) , 

A NPI11 > ,FACE<39,65) 

COMMON/DISTRN/CMMINI 10 ) , CMMAXI 10),CMAVG( 10 ) , CDPOP 1 ( 10 ) , CDP3P2 ( 10 ) , 
A CDPQP31 10) ,THTMIN( 10 ) , THTPLS 1 1 0 ) ,RMM IN ( 8 ) , RMMAX ( 8 ) »RMAVG( 8) , 

B RDPOPK 8 ) , RDPOP2 ( 8 ) * RDP0P3I 8 ) , TMM f N . TMMAX , TMA VG , SMAVG , TDP1 , 

C TDP2»TDP3*QLOC(10,20) 

COMMON/READTR/STATI, STATO,LCNTFL,LDPAR 
C 

C COMMON BLOCK FOR TRANSFERRING SPECIFIC DISTORTION 
C PARAMETERS TO THE OUTPJT SUBROUTINE 'LNPOUT* 

C 

COMMON/PARTRN/PSPECI 20) 

DIMENSION QD< 10), DTHETAI 10,20) 

DIMENSION RIDC(5),RIDR(5) 

LOGICAL STATI , S T ATO, L CNT FL , L DP AR 
C 

C CALCULATION OF KD2 
C 

N1 = 1 
N2 = NR 

I F ( STATI ) Nl=2 
IF(STATO) N2=NR-1 
D I NVS= 0. 0 
DO 26 I =N1 , N2 

DINVS=DINVSM 1.0/RADL0C( I) ) 

26 CONTINUE 

DO 27 I =N1 , N2 . 

27 QD ( I ) = ( 1 . O/RADLOC ( I ) ) /DI NVS 
RK02 = 0 . 0 

DO 28 I =N 1 » N2 

28 RKD2 =RKD2+ ( CDPOP3( I ) *THT MI N( I ) *QD( I 1*100.0) 

C • 

C CALCULATION OF K-THETA 
C 

DO 36 I =N1 , N2 
J 1 = NP ( I ) 

DO 36 J=1,J1 
K1=J-1 

IF(Kl.EQ.O) K1=J1 
Al=ANGLOC ( I , K1 ) 

I F ( { J-D.EQ.O) Al = Al-360.0 
K2= J+l 

IFIK2.GT.J1) K2= 1 
A2 = ANGLOC ( I ,K2) 

I F ( ( J4-D.GT.J1) A2 = A 2 + 360. 0 
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DTHETAI I, J)=(A2-Al)/2.0 
36 CONTINUE 

RGAM=. 285714 

QAV=(SMAVG/RGAM)*< ( ( ABS ( TM AVG/ SMAVG ) > **RGAM ) - 1 . 0 ) 

RKTHTO.3 
DO 35 I =N l » N2 
RISIN=0.0 
RICOSO.O 
J 1 = NP ( I ) 

DO 34 J=1,J1 

PH I = ANSLOC ( I , J I *.01745 3 

RISIN=RISIN-*-( (P( I»J) /CMAVG ( I))*SIN(PHI) *DTHETA( I , J) ) 
RICOS=RICOS+((P(I , J) / CMAVG ( I > )*COS{ PHI ) *DTHETA ( I * J > ) 

34 CONTINUE 

RKTHT=RKTHT+QD( I ) *( S QRT ( ABS( ( (ABS(RISIN) ) **2 .0 ) +t ( ABS( RICOS) ) 

A **2.5>>)) 

35 CONTINUE 

IFtQAV.EQ.'O.O) QAV=1.D 
RKTHT=(RKTHT/(180.‘0*QAV))*TMAVG 

CALCULATION OF IDC-MAX AND IDR-MAX 

NCHK=NR 

IF(STATI) NCHK=NCHK-1 
IF(STATO) nchk=nchk-i 

I F ( NCHK. NE. 5 ) WRITE(6,100) * • ' 

IF< NCHK.NE.5) WRITE(6,101) 

IFINCHK.NE.5) GO TO 40 * 

Jl = l 

J2 = NR * 

I F ( STATI ) Jl=2 
IF(STATO) J2=NR-1 
DO 41 J=J1,J2 

RIDCI J)=(CMAVG< J)-CMMIN( J) )/TMAVG 
41 RIDR(J)=(TMAVG-CMAVG< J) I/TMAVG 
RIDCXl=(RIDC(Jl)+RIDC(Jl+l))/2»0 
RIDCX2=(RIDC( Jl+3)+RIDC< Jl+4) 1/2.0 
RIDCMX=AMAXi ( RIDCX1, RIOCX2 ) 

RIDRMX=AMAX1 (RIDRIJ1 ) * R I DR l Jl+4) ) 

40 CONTINUE 

SPECIFIC DISTORTION PARAMETERS STORED IN ARRAY »PSPEC' 

IN PREPARATION TO BEING TRANSFERRED TO 'LNP0UT' 

PSPECI 1 )=RKD2 
PSPEC ( 2 ) =RKTHT 
PSPECI 3)=RIDCMX 
PSPECI 4)=RIDRMX 

100 FORMAT (10X*4DHI DC-MAX AND IDR-MAX CANNOT BE CALCULATED) 

101 FORMAT (10X»51HTHERE ARE NOT EXACTLY 5 RINGS OF TOTAL MEASUREMENTS) 
RETURN 

END 
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Subroutine PFIX 


When PFIX(. TRUE. ) is called by CFACE, QLOC (see appendix D) is calculated. 

When PFIX(. FALSE. ) is called by CFACE, PFIX transforms the array of input measure- 
ments (P) into an array of the proper mapping parameter (controlled by the value of 
WHICH) and writes it out. Also the call PFIX(. FALSE. ) by CFACE allows PFIX to step 
the input measurements outward by one ring if LCNTFL= .TRUE. . 

"NO STATIC RINGS, RESULTS HAVE NO MEANING" is a message that is written 
in two instances: First, if PFIX(.TRUE.) has been called by CFACE, the program will 
write the message and continue even though it is unable to calculate QLOC. If the user 
ever decides to modify the program to calculate QLOC without static -pres sure inputs, be 
sure it is done in the subroutine DISPAR. Second, where WHICH has been set to 3 
(WHICH=3, Mach no.) and there are no static-pressure input measurements. In this in- 
stance program execution will be terminated. 


SUBROUTINE PFIX(LQ) 

SUBROUTINE TO TRANSFORM INPUT MEASUREMENTS 
INTO MAPPING PARAMETERS 

SUBROUTINE ALSO CONTROLS SOME OUTPUT 

COMMON/READTR/STATI, STATO.LCNTFL.LOPAR 
COMMON/WRANGE/XBAR,RMAV, WHICH, RANGE, LPCENT,RMACH,DELX 

NP, RADLOC, ANGLOC, ANO P HAVE THE FIRST 

ARRAY DIMENSION EQUAL TO 11. 10 OF THE 

11 POSITIONS ARE AVAILABLE FOR USF. THE 

REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 

C0MM0N/DATTRN/P(11,20),NR,RA0L0C< li >, ANGLOC! 11,20) , 

A NP(ll) ,FACE(39,65) 

COMMON/DISTRN/CMMINI 10),CMMAX< 10),CWAVG( 10),CDPOP1 ( 10) ,CDP0P2( 13) , 
A C0P0P3 ( 10 ) , THTMI N ( 10 ),THTPLS< 10) ,RMMIN( 8 ) , RMMAXI 8 ) , RMAVG ( 8 ) , 

B RDPOPK 8 ) ,RDP0P2 ( 8 ) , RDP0P3 ( 8 ) , TMM I N, TMM AX , TMA VG, SM A VG , TDP1 , 

C TDP2,TDP3,QLOC( 10,20 ) 

COMMON/FINI/LQUIT 
DIMENSION NRFIT(2),Z(2),F(2) 

DIMENSION P S A ( 2 ) 

INTEGER WHICH 

LOGICAL LPCENT , LQU I T 

LOGICAL STATI,STATO,LCNTFL,LDPAR 

LOGICAL LQ 

IF(.NOT.LQ) GO TO 71 

DO 70 1=1,10 

DO 70 J=1 , 20 

QLOC ( I , J ) =0 . 0 

70 CONTINUE 

71 CONTINJE 
IF(LQ) GO TO 55 

SECTION THAT ENABLES I NT ERPOL AT I ON THROUGH THE 
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CENTER OF THE CONTOUR MAP 

I F ( . NOT . LCNT FL ) GO TO 50 
DO 51 1=1, NR 
NR1 =NR+2- I 

RADL0C(NR1 ) = P.ADL0C< NR1-1 ) 

K=NP ( NR1-1 ) 

00 52 J = 1 » K 

ANGLOC ( NR 1, J ) =ANG10C < NRl-1, J ) 

P(NR1, J)=P( NRl-1, J) 

52 CONTI N JE 
NP(NR1 )=NP(NR1-1 I 

51 CONTINUE 
N R = N R + 1 
RADLOC(1)=.01 
NP ( L ) = 3 
P S U M = 3 . 3 
K = NP ( 2 ) 

DO 54 1=1, K 
11=1-1 

IF( I 1.EQ.0) 1 1 = K 
Al=ANGLOC (2,11) 

IF ( ( I-D.EQ.O) Al = Al-360.0 

12 = 1+1 . : ' 

I F ( I2.ST.K) 12 = 1 - ' 

A2= ANGLOC (2,12) 

1 F ( (I+D.GT.K) A2 = A2 +360 .0 
DELTHTM (A2-A1)/ 2.0)/ 360.0 
PSUM=PSUM+(DELTHT*P( 2, I )) 

54 CONTINUE 

DO 53 1=1,8 
TH = ( I - 1 ) *45 
ANGLOC ( 1 , I ) =TH 
P(1,I)=PSUM 

53 CONTINUE 
50 CONTINUE 

CALCULATION OF MAXIMUM, MINIMUM, AND AVERAGE 
OF MEASUREMENTS . 

T MMAX= -500000.0 . 

T MM I N= +5000 00 * 0, J • 

PS=0.0 
RN=U • 3 
11 = 1 
I 2 = NR 

IF (STATU 11=2 
IF(LCNTFL) 11=2 
IF(STATO) I 2 =NR- 1 
DO 10 1=11,12 
J 1 = N P ( I ) 

DO 10 U=1,U1 
RN= RN+ 1.0 
PS=PS+P( I , J ) 

I F ( P( I ,J) .GT.TMMAX) TMMAX=P(I,J) 

I F { PCI ,U) .LT.TMMIN) TMM I N=P( I , J ) 

10 CONTINUE 
T M AVG=PS/RN 
C 
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C GO TO 11 FOR MACH NO. CONTOUR MAP 
C 

I F ( WHI CH. EO. 3 ) GO TO 11 

GO TO 14 FOR SPECIAL DISTRIBUTION CONTOUR MAP 
X-BAR */- N*( DELTA-X ) 

IFIWHICH.EQ.5.AND..N0T.LPCENT) GO TO 14 
MAPPING PARAMETERS GENERATED FROM INPUT MEASUREMENTS 
DO 12 1=1, NR 

J 1 =NP ( I ) - 

DO 12 J = 1 , J1 •' " ' ’ 

IF(WHICH.NE.l) P( I,J )=< ( (P< I, JI/TMAVG)-!. 01*100.0) 
IFIWHICH.EQ.l) P(I, J )=( (P( I, JJ/TMM AX) -1.0) *100.0 

12 CONTINUE 

SECTION TO WRITE OUT MAPPING PARAMETERS 

DO 19 J=1 , NR 
K=NP ( J ) 

IF( J.EQ.l. AND.LCNTFL ) GO TO 19 

IFIWHICH.NE.5) WRITE<6,100) RADLOCI J ) , < P ( J, I > , 1=1, K ) 
IFIWHICH.EQ.5) WRITE(6,lul) RADLOC IJ) , I P ( I , J ) , I =1 , K > 

19 CONTINUE 

MAPPING PARAMETERS PUT INTO PROPER FORM FOR SYMBOL SUBROUTINE 

IF(WHICH.EQ.l) GO TO 15 
DO 13 1=1, NR 
Jl^NPII) 

DO 13 J=1,J1 

PI I ,J> = (P< I , J)/( RANGE/20.3) J-10.0 

13 CONTINUE 
RETURN 

15 CONTINJE 

DO 20 1=1, NR 
J1=NP( I > 

DO 20 J=1,J1 

P(I ,J)=P(I,J)/( RANGE/ 20 . 0 ) 

20 CONTINJE 
RETURN 

MACH NO. SECTION 
11 CONTINJE 

ENTRY POINT FOR CALCULATION OF QLOC 
55 CONTINJE 

TEST FOR EXISTENCE OF STATIC RINGS 

I F ( .NOT.STATI . AND..NOT.STATO) GO TO 60 
GO TO 51 
60 CONTINJE 

WRI TE ( 5, 103 ) 

IFI.NOT.LQ) LQU I T= . T RUE. 
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RETURN 
61 CONTINUE 

INTERPOLATION TO FINO A STATIC PRESSURE FOR 
EACH TOTAL PRESSURE MEASUREMENT 

STATIC PRESSURE WILL BE A RADIUS WEIGHTED AVERAGE 
OF STATIC PRESSURE ON INNERMOST RING AND STATIC PRESSURE 
ON OUTERMOST RING, AT ANGLE THETA (IF EITHER RING IS 
MISSING, THEN IT WILL BE EXCLUDED FROM THE AVERAGE) 


I 1 = 1 
I 2 = NR 

I F ( ST AT I ) 11 = 2 
IF(STATO) I 2 = NR-1 
RMSUMO.O 
, KMSUMO 
DO 30 1=11,12 
J1=NP( I ) 

DO 30 J=1,J1 
PSA(1)=0.0 
PS A ( 2 ) =0. 0 
KCT = 0 
NR1=NR-1 

DO 31 K = 1 * NR, NR1 
KCT = KCT + 1 

I F ( .NOT.STATI.AND.K.EQ.l ) GO TO 31 
I F ( .NOT.STATO.AND.K.EO.NR) GO TO 31 
LI = NP ( K ) 

IF(Ll.NE.l) GO TO 37 
PSA(KCT)=P(K,1) 

GO TO 31 
37 CONTINUE 

I F ( ANGLOC ( I ,J).GT ,ANGL0C(K,L1) ) L=1 
IF(ANGL0C(I,J).GT.ANGL0C(K,Ll) ) GO TO 33 
DO 32 L=1 , L 1 

IF(ANGL OC(K,L).GE. ANGLOC (I, J)) GO TO 33 

32 CONTINUE 

33 CONTINJE 
NRF I T ( l ) =L- 1 
NRF I T( 2 ) =L 
DO 34 M= 1 , 2 

IF(NRFIT(M).LT.l) GO TO 35 
IF(NRFIT(M) .GT.L'l) GO TO 36 
N=NRF I T ( M ) 

Z ( M ) = ANGLOC ( K, N ) 

F(M)=P(K,N) 

GO TO 34 

35 CONTINJE 
N=NRF I T ( M ) +L1 
Z(M)=ANGLOC(K,N)-360.O 
F(M)=P(K,N) 

GO TO 34 

36 CONTINJE 

N = NRFI T ( M ) — L 1 

Z(M) =3S0.0*ANGL0C(K,N) 

F (M)=P(K,N) 

34 CONTINJE 

Z0= ANGLOC ( I , J) 

PSA ( KCT ) = 
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A f(( F f 2 ) -F ( 1 ) )/(Z( 2 ) -Z (1) ) )*IZO-Z( 1) > )*F( 1 ) 

31 CONTINJE 

RADIUS WEIGHTED AVERAGE STATIC PRESSURE 
RADIUS=RADLOC( I ) 

IF (PSA(l) .EQ.0.0 ) RADIUS =RADLOC(NR ) 

I F ( PSA (2) .EQ.0.0) RADIUS =RADLOC( 1 > 

PS=PSA(l)+( ( RADIUS-R ADLOC< 1) ) / ( RADLDC ( NR ) -R ADLOC ( 1) ) )* 

A (PSA(2I-PSA(1) ) 

MACH NO. IS CALCULATED FOR EACH INPUT MEASUREMENT LOCATION 
G AM= 1 . 4 

IF(PS.EO.O.O) PS=0. 00001 
IF(.NOT.LQ) GO TO 59 

QLOCII ,J) =( GAM/2.0 )*PS*SQRT( ABS(< 2. D/( GAM-1.0) )* 

A ( ( ABS ( P ( I, J) )**( (GAM-1.0) /GAM) >-1.0) ) ) 

59 CONTINJE 

IF(LQ) GO TO 63 
PTEMP=?(I,J) 

P ( I , J)=SQRT(ABS( (2.0/ (GAM-1.0) ) > 

A *({ (ABS(P(I»J)/PS) ) **( ( GAM-1.0) /GAM) )-1.0) ) ) 

IF(PTEMP.LE.PS) P(I,J)=0.0 
RMSUMsRMSUM+P ( 1 1 J ) 

KMSUM=<MSUM+1 
63 CONTINJE 
30 CONTINJE 

IF(LQ) RETURN 

CALCULATION OF AVERAGE MACH NO. 

RMA V=RMSUM/ FLOAT (KMSUM) 

SETTING MACH NO. = 0 ON RINGS OF STATIC PRESSURE 

I F ( . NOT . STAT I ) GO TO 38 
1 1 = MP ( 1 ) 

DO 40 1=1,11 

40 P(1 , I ) =0.0 

38 CONTINJE 

I F ( . NOT.STATO) GO TO 39 
I 1=NP( NR) 

DO 41 1=1,11 

41 P( NR, I )=0.0 

39 CONTINJE 
DO 43 1=1, NR 
J1=NP( I ) 

I F ( I .EQ.l. ANO.LCNTFL ) GO TO 43 
WRITE (5,102) RADLOC( I),(P(I,J),J=1,J1) 

43 CONTINJE 
DO 44 1=1, NR 
J1=NP( I ) 

DO 44 J=1,J1 

P ( I , J ) = ( P ( I , J ) /RMACH ) *( -20 .0 ) 

44 CONTINJE 
RETURN 

C 

C SECTION TO CALCULATE MAPPING PARAMETER FOR MAP 
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C OF X-BAR */- N*!L)ELTA-X) 

C 

14 CONTINUE 
C 

C CALCULATION UF X-RAR 
C 

XBAR=FLOAT(JFIX((( TMM AX + TMM I N ) / I 2.0*DELX ) ) +0.5) )*DELX 
DO 18 1=1, NR 
J1=NP( I) 

DO 18 J=1,J1 
C 

C PUTTING MAPPING PARAMETER INTO PROPER FORM FOR SYMBOL ROUTINE 
C 

P(I » J ) = ( (P( I »J)-XBAR)/DELX)-10.3 
18 CONTINUE - 

RETURN 

100 FORMAT! 5X,F 5.2, 4X, ID ( 2X,F8.4)/14X,13(2X,F8.4)) 

101 F0RMAT(5X,F5.2,4X,5( 5X, F8. 4) /3( 14X, 5! 5X, F8.4)/ ) ) 

102 FORMAT (5X,F 5.2, 4X, ID (4X , F6 . 4 ) / 14X , ID ! 4X , F6. 4 ) } 

103 FORMAT < IX, 4JHN0 STATIC RINGS, RESULTS HAVE NO MEANING) 

END 


Subroutine MAINLP 


MAINLP is the driving subroutine for the contour map interpolation. MAINLP 
calls INTERP and SYMBLE, in that order, for every r, 9 location in the contour map 
output grid. 


SUBROUTINE MAINLP 
C 

C DRIVING SUBROUTINE TO INDEX THROUGH THE GRID POINTS 
C IN THE INPUT MEASUREMENT PLANE 
C 
C 

C NP, RADLOC, ANGLUC, AND P HAVE THE FIRST 
C ARRAY DIMENSION EQUAL TO 11. ID OF THE 
C 11 POSITIONS ARE AVAILABLE FOR USE. THE 
C REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 
C 

COMMON/DAT.TKN/P! 11,20 ) ,NR, RADLOC! 11 ) , ANGLOC! 11,20), 

A NP(ll) ,FACE(39, 65) 

COMMON/ SC ALE/F ACTOR 
COMMON/ I NTRN/R, THETA, PV ALU E' 

DATA BLANK/ 1 H /,PLUS/1H+/ 

DO 10 10=1,39 
DO 1:0 J0 = 1 , 65 
C 

C CALCULATION OF R AND THETA 
C 

Y=( FLOAT (2D-IU) )*FACTOR 
X = ( ! FLOAT! JO-33) )/ 1 . 66 5 ) *F ACTOR 

r=sqrt(x*x+y*y) 

IF! X.EQ.O.D) GO TO 26 
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OMEGA=ATAN< ABS(Y/X) ) 

GO TO 27 

26 CONTINUE 
OMEGA= 1.5708 

27 CONTINUE 
C 

C CHANGING FROM CONVENTIONAL R-OMEGA SYSTEM TO THE 
C R-THETA SYSTEM USED IN THE INTERPOLATION PROCESS 
C 

IF(X.GE.O.O.AND.Y.GE.O.O ) THETA=1 . 5708-OMEGA 
IF(X.GE.O.D.AND.Y.LT.O.O > THETA=1 . 5708+0MEGA 
IF(X.LT.O.O.AND.Y.LT.O.O) THFTA=4. 71239-OMEGA 
IFIX.LT.O.O.AND.Y.GE.O.O ) THETA=4.71239+0MEGA 
I F( THETA. GE. 6. 28318) THETA=6. 28318 
I F { THETA.LT. 0.0) THET A = 0.0 
THETA=THETA*(360.0/6.28318) 

IFIR.LT. RAD LOCI l).OR.R.GT.RADLOC(NR) ) GO TO 9 
C 

C CALLING INTERPOLATION AND SYMBOL SUBROUTINES 
C 

CALL INTERP 
CALL SYMBLE 
FACE ( I 0 » JO ) =PVALUE 
GO TO 10 
9 CONTINUE 

FACE(IO,JO)= BLANK 
10 CONTINUE 

FACE(2D,33)=PLUS 

RETURN 

END 


Subroutine INTERP 


INTERP does the actual interpolation in the input measurement plane. It is fed 
values of r and 0, and it returns an interpolated value at that position. INTERP is 
first called by DISTRT when the simple radial distortion parameters are being cal- 
culated. At this time, interpolation is done using the input measurements in the array 
P. By the time INTERP is called by MAINLP, the array P contains mapping param- 
eters instead of input measurements, and the interpolation is done on the mapping 
parameters. 


SUBROUTINE INTERP 
C 

C INTERPOLATION SUBROUTINE 

C 

C 

C NP, RADLOC, ANGLOC, AND P HAVE THE FIRST 
C ARRAY DIMENSION EQUAL TO 11. 1C OF THE 
C 11 POSITIONS ARE AVAILABLE FOR USE. THE 
C REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION. 
C 

COMMON/DATTRN/PI 11,20), NR, RADLOC (11), ANGLOC! 11,20), 

A NP(ll) .FACE (39, 65-) 


o n o o o noon 


COMMON/ INTRN/R, T Hfc TA , RVALUE 

DIMENSION NFIT{2),NRFIT<2),RADFIT(2),Z(2),F(2) 

I F ( NR. LT. 2 ) WRI TE ( 6» 100 ) 

100 F0RMAT(1X,A9HY0U HAVE LESS THAN 2 RINGS, PROGRAM WILL NOT WORK) 
NR 1 = NR-1 

DETERMINE THE CORRECT TWO RINGS TO 
USE IN THE INTERPOLATION PROCESS 

DO II 1=1, NRI 
11 = 1 

IF(R.GE.RAOLOC( 1 1 > . ANO.R.LT. RADLOCI Il + l) ) GO TO 12 

11 CONTINJE 

12 CONTINUE 
NF I T ( 1 ) = 1 1 
NFIT(2) = IH-l 

SET-UP FOR CIRCUMFERENTIAL INTERPOLATION 
DETERMINE THE TWO MEASUREMENTS ON EACH RING TO USE 
IN THE CIRCUMFERENTIAL INTERPOLATION PROCESS 

DO 15 11=1,2 
K=NFIT< II ) 

J 1 = N P ( K ) 

IF(Jl.NE.l) GO TO 21 
RADFITI II)=P(K,1) 

GO TO 15 
21 CONTINJE 

IFITHETA.LT . ANGLOC ( K , 1) ) GO TO 19 

DO 16 J=1,J1 

J2=J 

I F ( J.EQ.J1) GO TO 17 

IF( THETA. GE.ANGLOC(K,J). AND.THETA.LT. ANGL0C(K,J+1 )) GO TO 17 

16 CONTINJE 

17 CONTINJE 
GO TO 20 

19 CONTINJE 
J2 = Ji 

20 CONTINUE 
NRF I T( 1 ) = J2 
NRF I T ( 2 ) = J2+ 1 
DO 18 1=1,2 

IF(THETA.LT.ANGL0C(K,1) ) GO TO 7 
I F ( NRF IT ( n.LT.l) GO TO 23 
IF(NRFIT( D.GT.Jl) GO TO 2 A 
L=NRFIT( I ) 

Z ( I ) =ANGLOO( K,L) 

FII ) = P ( K » L ) 

GO TO 18 ‘ 

7 CONTINJE 

IFLNRFITII I.LE.Jl) GO TO 6 
L=NRFIT(I)-J1 
Z ( I ) =ANGLOC ( K, L ) 

F ( I ) =P ( K, L ) 

GO TO 18 
6 CONTINJE 
L=NRFI T ( I ) 

Z( I )=ANGLOC(K,L)-360.0 
F { I) =P(K,L) 


28 



GO TO 18 

23 CONT IN JE 

L = NKF I T ( I )+Jl 
Z ( I >=ANGLOC(K, LJ-36D.0 
F ( I ) =P( K» L ) 

GO TO 18 

24 CONTINUE 

L = NRFIT( I )-Jl 
Z ( I )=3SO.O+ANGLOC(K,L) 

F(1)=P(K»L) 

18 CONTINUE 
C 

C LINEAR INTERPOLATION IN THE CIRCUMFERENTIAL DIRECTION 
C 

RAOFITI 1 1 ) = ( (F(2)-F( 1 ) )*( (THETA-ZI 1 ) )/( Z< 21-ZUH ))+F(l> 
15 CONTINUE 
C 

C LINEAR INTERPOLATION IN THE RADIAL DIRECTION 
C 

DO 30 1=1,2 
L=NF I T ( I ) 

Z ( I ) =RADLOC( L ) 

F ( I ) =RADFIT ( 1 ) 

30 CONTINUE 
C 

C FINAL INTERPOLATED VALUE 
C 

PVALUE = < ( F ( 2 I -F C 1 ) !*( (R-Z( 1) I/IZI2I-Z(1)I))+FI1I 

RETURN 

END 


Subroutine SYMBLE 


This subroutine determines the correct line-printer symbol to be placed on the 
contour map at each output grid position by the line-printer. 


SUBROUTINE SYMBLE 
C 

C SUBROUTINE TO SUPPLY THE SYMBOLS FOR 
C THE LINE PRINTER CONTOUR MAP 
C 

COMMON/ WRANGE/XBAR,RMAV, WHICH, RANGE , LPCENT , RMACH, DEL X 
COMMON/ I NTRN/R, THETA, BETA 
INTEGER WHICH 
LOGICAL LPCENT 

DIMENSION SYMBOL (23) , OUTSYMC 20 ) 

C 

C TWO LISTS OF SYMBOLS COVER ALL POSSIBLE TYPES OF CONTOUR MAPS 
C 

DATA SYMB0L/1H3,1H1, 1H2, 1H3, 1H4, 1H5, 1H6 , 1H7, 1H8, 1H9, 

A 1HA, IHB,1 HC, 1HD, 1HE, 1HF, 1HG, 1HH, 1HI , 1HJ/ 

DATA 0UTSYM/IH9, 1H8, 1H7, 1H6, 1H5.1H4, 1H3, 1H2, 1H1 , 1H3, 1H A , 1HB , 
A 1HC, 1HD, 1HE, 1HF,1HG, 1HH, 1HI , 1HJ/ 

IF(BETA.LE.-2U.3) BE T A.= - 19 . 9999 
IF( BETA. GE. 3.0 ) BETA=-0.0001 
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c 

C SECTION TO DETERMINE THE PROPER SYMBOL FOR THE CONTOUR MAP 
C 

DO 10 1=1,20 

UP=1-I 

DN=UP-1.0 

I F { BETA. LE, UP. AND. BETA. GT. ON) GO TO 11 

10 CONTINJE 

11 CONTINJE 

BET A = 0 JTS YM ( I ) 

I F ( WHICH. EQ. I.OR.WHI CH.EQ. 3) BETA= SYMBOL ( I ) 

RETURN 

END 


Subroutine LNPOUT 


LNPOUT controls the output of the contour map and the distortion parameter, in 
formation. In addition to the output of the simple distortion parameters, LNPOUT 
contains the user-written output formats for the user -calculated specific distortion 
parameters 


SUBROUTINE LNPOUT 
C ' 

C SUBROUTINE THAT CONTROLS THE OUTPUT 
C SOME OUTPUT CONTROLLED BY •CARDRD* AND *PFIX* 

C 

C - 

C NP, RADLOC, ANGLOC, AND P HAVE THE FIRST 
C ARRAY DIMENSION EQUAL TO 11. 1C OF THE 
C 11 POSITIONS ARE AVAILABLE FOR USE. THE 
C REMAINING POSITION IS NEEDED FOR INTERNAL MANIPULATION; 

C 

COMMON/DATTRN/P( 11,20 ) , NR, RADLOC ( 11 >, ANGLOC ( 11 ,20) , 

A NP( 111, FACE! 39, 65) 

COMMON/READTR/STATI, S T ATO, LCNTFL , L DP AR 
COMMON/WRANGE/XBAR,RMAV, WHICH, RANGE, LPC ENT, RMACH,DELX 
COMMON/DISTRN/CMMIN( 10), CMMA X ( ID) , CMAVG( 10) ,CDP0P1 ( 10 ) , CDPOP2 (10) 
A CDPOP3 ( 10) , THTMI N ( 10 ) , THTPL S ( 1 0 ) , RMM I N ( 8 ) , RMMAX ( 8 ) ,RMAVG<8), 

B RDPOPK8) ,RDP0P2 (8) ,RDPOP3( 8 ) , TMM I N, TMM AX, TM AVG , SMAVG , TOP 1 , 

C TDP2»TDP3,QLOC (10,20) 

COMMON/ HE ADG/T I TLE1( 80), TITLE2( 80) 

COMMON/PARTRN/PSPEC( 20) 

DIMENSION ROUT (19) 

INTEGER WHICH 
LOGICAL LPCENT 

LOGICAL STATI ,STATO,LCNTFL,LOPAR 
C 

C LOGIC TO DETERMINE WHICH KIND OF MAP 
C IS TO BE PRINTED 
C 

IF(WHICH.EQ.5) GO TO 11 

if(whi:h.eq.<,) go to 12 

IF(WHi:H.EQ.3) go to 13 
I F ( WHI CH. EQ. 2 ) GO TO 14 
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OUTPUT SE:TI ON FOR ( PI-PMAX )/PMAX MAP 
00 21 1=1,19 

ROUT ( I ) =-FLUAT ( I )*(R4NGE/20.0) 

21 CONTINJE 
WRITE! 6*100) 

WRI TE ( 6 * 1 90 ) TITLE1 
WRITE(6,190) TITLE2 
WRI TE( 6,105 ) 

WRI TE( 6,106 ) 

WRITE! 6, 108) ( ROUT (I), 1=1, 19) 

WRITE! 6,109) - 

WRITE16.132) 

WRITE(6,138) T MM AX 
GO TO 15 

OUTPUT SECTION FOR ( Pl-PAVG )/PAVG MAP 

1A CONTINJE 

DO 22 1=1,19 

ROUT! I )=FLQ AT! ID-I )*( RANGE/ 20.0) 

22 CONTINJE 

WRI TE! 6,100 I 
WRI TE! 6,190 ) TITLE1 
WR I T E ( 6 , 1 90 ) TITLE2 
WRITE! 6,105) 

WRITE!6,1D7) 

WRI T E ( 6 , 1 10 ) (ROUT! I ), 1 = 1, 19) 

WRI TE! 6,109) 

WRITE! 6,128) 

WRITE! 6,137 ) TMAVG 
GO TO 15 
13 CONTINJE 

OUTPUT SECTION FOR MACH NUMBER MAP 
DO 20 1=1,19 

ROUT ( I )=(RMACH/20.0)*FL0AT! I ) 

20 CONTINUE 

WRI TE! 6,100 ) 

WRI T E ( 6 * 190 ) TITLE1 
WRITE(6,190) TITLE2 
WRITE! 6,105) 

WRI TE ( 5 * 1 06 » 

WRITE(6,111 ) ! ROUT (I ) , 1 = 1, 19 ) 
WRITE(6,133) 

WR ITE(6»1A1 ) RMAV 
GO TO 15 
12 CONTINUE 

OUTPUT SECTION FOR ( T I - T AVG ) / T AVG MAP 
DO 23 1 = 1 ,19 

ROUT! I ) =FLOAT ( 10-1 )*( RANGE/20. 0 ) 

23 CONTINJE 

WR I T E ( 6 » 1 00 ) 

WR I TE ! 6 , 1 90 ) TITLE1 
WRI TE! 6,190 ) TITLE2 
WRI T E ( 6 , 1 90 ) TITLE2 
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WRITEC 5,105) 

WRI TE( 5,107 ) 

WRI TE ( 5 , 1 1 2 ) C ROUTC I ) , 1=1, 19) 

WRI TEC 5,109 ) 

WRI TECS, 129) 

WRI TEC £>,139) TMAVG 
GO TO 15 

OUTPUT SECT l ON FO MAP WITH SPECIAL INPUT 
11 CONTINUE 

I F ( .NOT.LPCENT) GO TO 17 
OUTPUT SECTION FOR ( X I - X AVG )/ X AVG MAP 
DO 24 1=1,19 

ROUT ( I )=FLOAT< 10- I )*< RANGE/ 20. 0 ) 

24 CONTINUE 

WR I T E ( b , 1 00 ) 

WR I T E ( b , 1 90 ) TITLE1 
WRITE! b, 190) TITLE2 
WRITEC 5,105) 

WRI T E ( 5 , 1 07 ) 

WRI TE ( 5, 1 13 ) CROUTCI), 1=1,19) 

WRITEC 5,109) 

WRI TE C 5 , 1 31) 

WRITEC5.130) 

WRITEC 5, 140) TMAVG 
GO TO 15 
17 CONTINUE 

OUTPUT SECTION FOR MAP OF (X-BAR) +/- N(DEL-X) 

WRI TEC 5,100 ) 

WRITEC5, 190 ) TITLE1 
WRITE(5,190) TITLE2 
WRITE(5,105) 

WR I T E ( 5 » 1 07 ) 

WR I TE ( 5 , l 1 4 ) 

WRITEC b, 115) 

WRITEC 5,116) 

WRITEC 5,117 ) 

WRI TEC b, l 34 ) 

WRI TEC 5,135 ) 

WRI TEC 5,142) XBAR 
WRI T E ( 5 , 1 43 ) DELX 
15 CONTINUE 

WRI TEC 5, 120) 

DO 10 1=1,39 

IFCl.NE.20) WR ITEC6, 101) < FACE l I , J ) , J=1 , 65 ) 
IFCI.E0.20) WRITEC 6, 121) C FACE ( I , J ) , J=1 , 65 ) 
10 CONTINUE 

WRITEC5.122) 

OUTPUT SECTION FOR SIMPLE DISTORTION PARAMETERS 

IF ( . NOT. LDPAR ) GO TO 16 
WRI TEC 5,100 ) 

WRI TEC 5,190) TITLEl 
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c 

c 

c 


WRITERS, 200 ) 
WKI TE! 6,201 ) 
WRITE ( 5,2 02 ) 
WRI TE ( 5,203) 
WRITE (5,204) 
WRI TE ( 5,205 ) 
WRI TE( 5,206 ) 
WR I TE ( 5 » 207 ) 
WRI T E ( 5,203 ) 
WRITE(5,209) 
WRI TE( 5,210) 
WRITE! 5,211 ) 
WRITE! 5,212 ) 
WRI T E ( 5 , 2 1 3 ) 
WRITE! 6,214 ) 
WRITE! 5,215 ) 
WRI TE! 5,216) 
WRI TE(S,217) 
WRI TE! 5,218) 
WRI TE! 5,219 ) 
WRI TE! 5,220) 
WRI T E ( 5 , 2 2 1 ) 
WRI TE! 5,222 ) 


( THTMIN! I ) , 1=1, ID) 
( THTPLS ( I ) , 1=1 , 10 ) 
( CMMIN! I ), 1=1, 10) 

( CMMAX ( I ), 1=1, 10) 
(CMAVG! I), I =1,10) 

( CDPQP1 ( I >, 1=1, 10) 
( CUPOP2 ( I ), 1=1, 10) 
( CDP0P3 ( I ), 1=1, 10) 


(RMMIN! I ) , I =1 , 8 ) 

( RMMAX ( I ) , 1=1, 8 ) 

( RMAVG! I ) , I =1, 8 ) 
(RDP0P1! I ), 1=1,8) 
(RDP0P2! I ), 1=1,8) 
( RDP0P3 ( I), 1 = 1,8) 

TMMI N, T DPI 
T MMAX , T DP2 
T MAVG, T DP3 
SMAVG 


OUTPUT SECTION FOR USER-DEFINED DISTORTION PARAMETERS 


WR I TE ( 6 , 1 00 ) 

WRI TE ( 5,190) TITLE1 
WR I rE(5,190) TITLE2 
WRITE! 5,10?) 

WRI TE ( 5,300) 

WR I T E ( 6,301 ) PSPEC(l) 

WR I TE ( 6 , 302 ) PSPEC<2) 

WRI TE! 5,303 ) PSPEC(3) 

WRITE( 5,304) PSPEC ( 4 ) 

WRITE! 6, 100) 

16 CONTINJF 

100 FORMAT ( 1H1 ) 

101 FORMAT ( 1X»33X,65A1 ) 

102 FORMAT! 1H0) 

105 FORMAT ( IX/ / 5 IX , 22HKE Y TO MAPPING SYMBOLS/) 

106 F ORM AT ( 2 X , 6HB0RDER, 1 1 X , 3H0/ 1 , 3X , 3H1 / 2, 3X, 3H2/3 , 3X , 3H3/4, 3X , 

A 3H4/5,3X, 3H5/6.3X, 3H6/7, 3X,3H7/8, 3X, 3H8/9, 3X , 3H9/A , 3X , 3H A /B , 

B 3X,3HB/C,3X,3HC/D, 3X , 3HD/E , 3X , 3HE /F , 3X , 3HF /G , 3X , 3HG/H , 3X , 3HH/ 1 , 
C 3X.3HI/J/) 

107 FORMAT (2X ,5HBURDER,1 IX, 3H9/ 8,3X,3H8/7,3X,3H7/6,3X,3H6/5,3X,3H5/4, 
A 3 X , 334/ 3, 3X, 3H3/2, 3X, 3H2/ 1 , 3X, 3H1/0, 3X, 3H0/A, 3X, 3HA/B, 

B 3X,3HB/C,3X,3HC/D, 3X , 3HD/E , 3X , 3HE/F , 3X , 3HF/G , 3X , 3HG/H, 3X , 3HH/I , 
C 3X.3HI/J/ ) 

108 FORMAT (2X,14H(P I -PMAX)/PMAX, 20! IX, F 5.1) ) 

109 FORMAT ( 2 X » 1 OHI N PERCENT) 

110 FORMAT (2X,14H( PI -PAYG J/PAYG, 20! IX, F 5.1) ) 

111 FORMAT (?X,1 1HMACH N JMBER , 3X, 20 ( 2X , F 4. 3 ) ) 

112 FORMAT (2X» 14H( TI-TAYGI/T AVG* 20! 1X,F5.1) ) 

113 FORMAT ( 2X , 1 4H ( x I — X AY G ) /X AYG, 20 ( IX, F 5.1) ) 

1 1 A FORMAT ( 2X , l 4HMUL.T I PL I ER FOR, 4X , 2H*9, 4X , 2H* 8 , 4X , 2H*7 , 

A 4X,2H+6,4X,2H + 5,4X ,2H+A f 4X , 2H*3 , 4X , 2H+2 , 4X, 2H+ 1 ,5X,1H0, 

B 4X,2H-1,4X,2H-2,4X,2H-3,4X,2H-4,4X,2H-5,4X,2H-6,4X, 

C 2H-7,4X,2H-8,4X,2H-9) 
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115 FORMAT ( 2 X ,1 faHDE t T A-X O/E) 

116 FORMAT) 2X , 1 3HBORDER CORRESPONDS) 

117 FORMAT (2X,23HTU ( X-B AR ) + ( - 4*DEL T A-X ) ) 

120 FORMAT) 1X/66X,16H+ 360/O DEGREES/) 

121 FORMAT (1X»25X,6H270 +, 2X, 65 A1 , 2X, 5H+ 90) 

122 FORMAT (1X/66X,6H+- 180) 

128 FORMAT ( IX////51X»23H ( PI-PAVGI/PAVG PROFILE) 

129 FORMAT) 1 X// / /5 1 X , 23H ( T I -TAVG l/TAVG PROFILE) 

130 FORMAT ( 5 1 X, 2 3H ( XI-XAVO/XAVG PROFILE) 

.131 FORMAT) 1X///51X, I3HSPECIAL INPUT) 

132 FORMAT (1X////51X,23H) PI- PM AX ) / PMAX PROFILE) 

133 FORMAT) 1X////51X.16HMACH NO. PROFILE) 

134 FORMAT )51X,2DHSPECIAL DISTRIBUTION) 

135 F OR MAT )51X»28HX-BAR PLUS/MINUS NMOELTA-X)) 

137 FORMAT)51X, 17HAVERAGE PRESS JRE = , F7. 3 ) 

138 FORMAT ( 51 X, 17HMAXI MUM PR ES SURE= , F7 . 3 ) 

139 FORMAT )51X»20HAVERASE TEMPERATURE=, F7.2 ) 

140 FORMAT ( 51 X , 1 6H AVER AGE SPEC I AL = , F 1 1 „ 5 ) 

141 FORMATI51X..17HAVERAGE MACH N0.=,F6.4) 

142 FORMAT (51X,6HX-BAR=»F11.5) 

143 FORMAT)51X,8HDELTA-X=,F9.4) 

190 F0RMAT(6X,8UA1) 

200 F0RMAT(1X//45X»37HCIRCUMFERENTIAL DISTORTION PARAMETERS//) 

201 FORMATI3X,8HRING NO. »17X,1H1,7X,1H2,7X, 1H3,7X, 

A 1H4,7X,1H5,7X,1H6»7X,1H7»7X,1H8,7X,1H9,6X»2H10//) 

202 FORMAT(3X,6HTHTMIN,14X,10)2X,F6.2>/) 

203 F0RMAT(3X»6HTHTPLS,14X, 10) 2X,F6.2)/ ) 

204 FORMAT ( 3X,5HCMMIN»15X» 10 ( 1 X , F7.3)/ ) 

205 FORMAT(3X»5HCMMAX,15X,10(1X,F7.3)/) 

206 FORMAT(3X»5HCMAVG,l5X,iO(iX,F7.3)/) 

207 FORMAT <3X, 1 9H < CMMAX-CMM I N ) /CMM AX » IX, 10) 2X.F6.5)/) 

■ 208 FORMAT )3X,19H( CMM AX'- CMM I N)/CMAVG»lX»10(2X»F6o5)/) 

209 F0RMAT)3X,19H(CMAVG-CMMIN)/CMAVG, IX, 10) 2X,F6.5)/ ) 

210 F0RMAT(1X///45X,28HRADIAL DISTORTION PARAMETERS//) 

211 FORMAT )3X,20HTHETA POSITION ( DEG ) , 5X, 1HD, 6X , 2H45 , 6X , 2H90 » 

A 5X,3H135,5X,3H18D, 5X,3H225,5X, 3H270, 5X,3H315//> 

212 F0RMAT(3X,5HRMMIN»15X*8( 1X»F7.3)/) 

213 FORMAT )3X,5HRMMAX,15X,8( 1X,F7.3)/) 

214 F0RMAT(3X,5HRMAVG,15X,8) 1X.F7.3)/) 

215 FORMAT ( 3X, 1 9H ( RMMAX-RMM IN) /RMM AX, iX,8(2X,F6.5)/) 

216 F ORMAT ( 3'X , 1 9H ( RMMAX-RMM I N) / 'RMAVG, IX , 8') 2X ; , F 6o 5 ) / ) 

217 F0RMAT)3X,19H(RMAVG-RMMIN)/RMAVG, IX , 8) 2 X, F6 . 5 ) / ) 

218 FORMAT) 1X///28X, 63HPARAMETERS PERTAINING TO THE ENTIRE ARRAY OF IN 
APUT MEASUREMENTS/) 

219 F0RMAT(38X»6HTMMIN=,F7.3»8X,5HTDPl=»F6o5) 

220 FORMAT!38X»6HTMMAX=,F7.3,8X,5HTDP2=»F6o5) 

221. FORM AT ( 38 X, 6HTMAVG", F7 . 3 , 8X » 5HTDP 3= , F6. 5 ) 

222 FORMAT) 39 X, 6 HSMAVG= » F7. 3 ) 

SECTION FOR FORMATS FOR THE OUTPUT 
OF USER-DEFINED DI STORTI ON. PARAMETERS 

300 FORMAT (1X//6X,34HUSER-DEFINED DISTORTION PARAMETERS///) 

301 FORMAT ( 1X/1 IX, 3HKD2, 8X.F11.5 j 

302 FORMAT ! 1X/11X ,7HK— THET A, 4X-, F 1 1 o 5 ) 

303 FORMAT) IX/ 11 X,9H) I DC ) -MAX, 2X , FI 1 . 5 ) 

304 FORMAT) IX/ 11 X,9H( I DR ) -MAX, 2X.F11.5) 

RETURN 
END 
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APPENDIX B 


METHOD OF INTERPOLATION 

Given the array of input measurements as discussed in the section Instructions For 
Preparing Input, an interpolated value of the input measurement for any radial and cir- 
cumferential position (radius r and angle 0) in the measurement plane is desired. 
Figure 2 shows a hypothetical array of input measurements. The user's array may 
differ from the example array in the number of rings and in the number of measure- 
ments on each ring. The ring radii and the angular positions of the measurements may 
also be different. 



Figure 2. - Hypothetical array of input measurements. 


The interpolation process used by the program is shown in figure 3. Figure 3(a) 
shows the circumferential interpolation process, and figure 3(b) shows the radial inter- 
polation process. 

Interpolation is done in the circumferential direction first. An interpolated value 
of the input measurement is needed at circumferential position 0 for each of two rings: 
one on either side of radius r. In the case shown in figure 2, rings J and J+l are used. 
If the user's array of input measurements contains less than two rings, a message to 
that effect will be printed out and program execution will be terminated. 

To find an interpolated value on each ring at angle 0, two measurements on each 
ring are used: one on one side of angular position 0 and one on the other side. For 
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example, ring J would use measurements 1 and 2 while ring. J+l would use measure- 
ments 8 and 1. If there are less than two input measurements on any given ring, the 
whole ring is assumed to have a constant value equal to the sole input measurement on 
that ring. 

A linear interpolation is then performed on both rings using the two input measure- 
ments (and their corresponding angular positions) on each ring. Referring to figure 3(a), 
this would result in obtaining interpolated values A and B for rings J and J+l, re- 
spectively. 

Finally, a linear interpolation is performed in the radial direction to arrive at the 
final interpolated value of the input measurement. The radial linear interpolation uses 
the two interpolated values (and their respective radii) resulting from the. two circum- 
ferential linear interpolations. Referring to figure 3(b), interpolated values A and B 
are used to linearly interpolate for the final interpolated value at radius r and angle 9. 

Linear interpolation has been proven to be sufficiently accurate for generating 
contour maps. This determination is based on comparing the computer generated radial 
and circumferential profiles with those actually obtained from full-scale engine testing. 


-Interpolated values 



Ring J 

Angular position 

(a) Circumferential interpolation. 


c 

° E 
oj £ 

— i-. 

> s 

E 




/-Final 

/ interpolated 
/ value 


:/ 
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,-r (See fig. 2) 
t 


r J 


r J+l 


Radial position 
(b) Radial interpolation. 


Figure 3. - Method of interpolation. 
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APPENDIX C 
SYMBOLS 

CM input measurement located on given ring 

D diameter of given ring of measurements 

NP number of measurements on ring 

NR total number of rings in array of input measurements 

OD outer duct diameter 

P pressure, total or static, kPa (psi) 

1 2 

Q dynamic pressure (— pu ) 

2 

RM input measurement located at discrete 9 position 

T temperature, total or static, K (°R) 

r length of radius vector in plane of input measurement, cm (in.) 

u flow velocity in axial direction, m/sec (ft/sec) 

X special input variable, user -defined 

AX increment of special input variable represented by each mapping symbol 

X mean value of special input variable rounded off to nearest AX 

(AP) rmS true root -mean-square value of time-variant pressure 

y ratio of specific heats 

9 reference angle in plane of measurement, measured clockwise from the 

12 o'clock position, deg 

0~ circumferential extent of largest region wherein the interpolated value of 

the measurement is below average value of input measurements on that 
ring, deg 

0 + defined as 360° -£®k’ where £ 0^. is summation of circumferential ex- 

tents of all areas on a given ring, wherein the interpolated value of the 
measurement is below average value of the input measurements on that ring 

p air density, kg/m^ (lbm/ft^) 

Subscripts: 

av average value of measurements 
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over whole array of input measurements 
denotes particular measurement on given ring 
denotes local value of measurement 
denotes particular ring 
maximum value of measurements 
minimum value of measurements 
static (pressure or temperature) 
total (pressure or temperature) 



APPENDIX D 


CALCULATION OF DISTORTION PARAMETERS 

As mentioned in the section OUTPUT, when the logical variable LDPAR is set to 
.TRUE. , simple circumferential and radial distortion parameters as well as specific 
distortion parameters are calculated. 

For each ring, the circumferential parameters consist of 

(1) Theta-minus (0"), the circumferential extent (in deg) of the largest area wherein 

the interpolated value of the measurement is below the average value of the 
input measurements on that ring (This parameter is in an array named 
THTMIN . ) 

(2) Theta-plus (B 4 ), defined as 360° - £9^ where £0^ is the summation of the 

circumferential extents of all of the areas on a given ring wherein the inter- 
polated value of the measurements is below the average value of the input 
measurements on that ring. (This parameter is in an array named THTPLS.) 

(3) The minimum input measurement value on a ring, CM min , stored in an array 

named CMMIN 

(4) The maximum input measurement value on a ring, CM,,,,,, stored in an array 

max 

named CMMAX 

(5) The average input measurement value on a ring, CM av , stored in an array 

named CMAVG 

(6) ( CM max " CM min^ CM max’ stored in arra y named CDPOP1 

(7) (CM max - CM min )/CM av , stored in array named CDPOP2 

(8) (CM„ V - CM._ )/CM„„, stored in array named CDPOP3. 

dv mm dv 

These parameters are calculated from input values. No interpolated values are used in 
calculating distortion parameters in the circumferential direction except in the case of 
9" and 0 + . 

At each of the eight discrete 0 positions described in appendix F (0°, 45°, 90°, 
etc.), the following simple radial distortion parameters are calculating using the cir- 
cumferential interpolations from each ring: 

(1) The minimum interpolated value at angular position 9, RM m ^ n , stored in an 

array named RMMEN 

(2) The maximum interpolated value at angular position 0, RM_ . stored in an 

nidX 

array named RMMAX 

(3) The average interpolated value at angular position 0, RM ay , stored in an array 

named RMAVG 

(4) (RM„,„„ - RM )/RM m , v , at angular position 0, stored in an array named 

max ill ill max 

RDPOP1 
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(5) (RM - RM • )/RM , at angular position 8, stored in an array named 

IXldA • illiil d.V 

RDPOP2 

(6) (RM„ , - RM . J/RM„. at angular position 0, stored in an array named 

a.v min 3,v 

RDPOP3 . 

Eight other parameters are calculated from the whole array of input measurements 
and are as follows: 

(1) The minimum total input measurement, TM called TMMEN 

(2) The maximum total input measurement, TM max , called TMMAX 

(3) The average total input measurement, TM ay , called TMAVG 

(4) The average static input measurement, SM„„, called SMAVG 

( 5 > (™max - ™min)/™max’ called TDP1 

(6) (™max ' ™min )/TM av’ called TDP2 

(7) (TM av - TM min )/TM av , called TDP3 

(8) Local values of dynamic head, stored in an array named QLOC. QLOC is 

dimensioned as QLOC(10, 20). (There is a value of for every correspond- 
ing input measurement. QLOC is not printed out but is available in subroutine 
DISPAR for use in calculating specific distortion parameters. Values of QLOC 
will only be calculated when the logical variable LDPAR is set equal to 
.TRUE. . Inner ring and/or outer ring static measurements are necessary, 
and of course the input measurements need to be pressures. If there are no 
static measurements, the Q/s will be set equal to 0.0.) 

The User may want to calculate a more complex, specific distortion parameter 
based on the array of input measurements. More often than not, a specific parameter 
will use at least some of these simple distortion parameters along with the basic 
geometry of the input array. Four examples follow showing how to calculate specific 
distortion parameters at an engine inlet. 

The user should place program steps only in the subroutine DISPAR and should be 
careful to not alter the input data in the common block PARTRN. Altering the program 
in any other subroutine could cause problems with the remainder of the program 
(especially the generation of the contour map) . ‘ 

The first example involves a distortion parameter introduced by Pratt & Whitney 
(see ref. 3) called KDg, which is defined by the following equation: 


KD 2 = 


NR 


J 


NR 


Y OD 
tl D J 
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where 


J particular ring of total pressure measurements 

0j circumferential extent, of largest single total -pressure depression below 

(P, „ ) T for a particular ring, deg 

L, ciV J 

(AP/P)j [j(P t ay - P t) min )/P t> av Jj for a particular ring, percent 
(P t ay )j average measured total pressure for particular ring 
(P t min )j minimum measured total pressure for particular ring 
Dj diameter of particular ring 

OD outer duct diameter (Note that OD is a constant and, since it appears in 

both the numerator and the denominator, it can be cancelled out.) 

NR number of rings of total pressure measurements (Remember that this is 

part of the input data.) 

Some of the quantities needed to compute KDg have already been calculated by the 
program. The value of 9“ for each ring is stored in THTMIN(J); (AP/P)j can be 
easily generated as CDPOP3(J)*100. 0; and Dj can be calculated simply as 
RADLOC(J)*2.0. The parameter KD 2 can then be generated by the following 
FORTRAN IV statements: 

SUM1-0.0 

SUM2=0.0 

Jl=l 

J2=NR 

IF (ST ATI) Jl=2 
IF(STATO) J2=NR-1 
DO 10 J=J1, J2 

• ' » 

SUM1=SUM1+1 . 0/(2 . 0* RADLOC(J)) 

10 SUM2=SUM2+THTMIN(J)*CDPOP3(J)*100. 0*(r. 0/(2.0*RADLOC(J))) 
RKD2=SUM2/SUM1 | 

The tests on STATI and STATO are necessary because no static -pressure rings are 
used in the calculation of KD 2 - j 

The second distortion parameter to be discussed was also introduced by Pratt & 
Whitney (see ref. 4) . The parameter is called Kg and is defined by the following 
equation: 
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1 


K e = 


1’ *.v 

Qay V .1 

“ D J 


The terms appearing in the equation for Kg are defined as follows: 
particular ring of total pressure measurements 

J±\ 1 maximum (AjVL^)j for L = 1, 2, 3, ... 


12 

dynamic head at plane of measurement pu ), calculated from average 

• * w 

total and average static pressures 
average of all total pressure measurements in input plane 
diameter of a particular ring 
number of rings of total pressure measurements 


t, av, f 


< A L>J = V a L + b L 


The Fourier coefficients a^ and b^ are defined as 


a L = 


i r p t < 9 > 
J P t, av, 

" rr 


cos L0 d0 all values from ring J 


i r p t^ 

b T = - / — si 

L 11 I P 

J. t,av,f 


sin L0 d0 all values from ring J 
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where 


P t (0) local total pressure measurement 

0 reference angle in measurement plane 

For most commonly encountered distortions, (A^/L^) = A^, and K ( 

written as 



The expression (A., ) can be written as 

1 J 


< a i> - 


where 


NP 


1 1802l-< P* 


Pf(«r) 


jTj t, av, f 


(cos 0j) A0j all values from ring J 


NP 


1 V P t<V 

b., = > (sin 0 T ) A0 T all values from ring J 

1 180/— i P f aw . 1 11 

j-j t, av, f 


and 


NP number of total pressure measurements on ring J 
I denotes particular measurement on ring J 

A0j l/2(0j + i - for particular ring J and 0j (A0j is in degrees) 

Hence, Kq can now be written as 


can be 
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The local total pressures are in the input measurement array P. Measurement I 
on ring J. is P(J, 1) . The value of P, f is already calculated and is called TMAVG. 

L y d V y X 

Again, Dj = 2. (j*RADLOC(J). The 0j's are stored in ANGLOC, and 0j for ring J is 
ANGLOC(J,I). The value of (P, „„) T is equal to CMAVG(J), and A0 T is equal to 
(ANGLOC(J,I+1)-ANGLOC(J,I-1))/2.0 for ring J. The value of NP simply equals NP(J) 
while Q„„ can be calculated from the equation 

dV 



where P s av j is the average of all static pressure input measurements and y is the 

ratio of specific heats. The value of P„ t is already calculated and is called 

s, <tv, I 

SMAVG. Note that, if there are no rings of static pressures, Kg cannot be calculated 
unless another method is used to arrive at a value of Q ay . One such method would be 
to input a reasonable value of average static pressure based on corrected airflow or . 
some other knowledge of the flow conditions at the input measurement plane. 

Now Kq can be calculated by the following FORTRAN IV program steps: 
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c 


c 


c 


calculation of qavg 

GAM=1.4 } a value of y was assumed 
RGAM=(GAM-1 . 0)/GAM 

QAVG =(SMAVG/RGAM)* (((TMAVG/SMAVG)** RGAM) -1 . 0) 
■CALCULATION OF DELTA -THETA 


Jl=l • 


J2=NR 


IF (ST ATI) Jl=2 
IF(STATO) J2=NR-1 


DO 10 J^Jl, J2 
NPJ=NP(J) 

DO 10 1=1, NPJ 


The summation from J=1 to J=N in the equation for Kq is 
done for the first total -pressure ring through the last total 
pressure ring. 


K1=I-1 

IF(Kl.EQ.O) K1=NPJ 
Al=ANGLOC(J, Kl) 
IF((I-1).EQ.0) A1=A1-360.0 


K2=I+1 


IF(K2.GT.NPJ) K2=l 
A2=ANGLOC(J, K2) 
IF((I+1).GT.NPJ) A2=A2+360.0 
10 DTHETA(J, I)=(A2-Al)/2.0 
— CALCULATION OF K -THETA 
SUM1=0.0 


SUM2=0 . 0 


DO 11 J=J1, J2 
RSIN=0.0 
RCOS=0 . 0 
NPJ=NP(J) 

DO 12 J=l, NPJ 

THETA=ANGLOC(J, 1)* 0.017453 

RCOS=RCOS+ (P(J, I)*COS(THETA)*DTHETA(j, I)) 

12 RSIN=RSIN+ (P(J,I)*SIN(THETA)*DTHETA(J,I)) 

SUM2=SUM2+(1 . 0/(RADLOC(J)*2. 0)*SQRT ((RCOS** 2 . 0)+(RSIN** 2 . 0)) 
11 SUM1=SUM1+(1 . 0/ (RADLOC(J)* 2.0)) 

RKTHT=SUM2/ (180. 0*QAVG*SUM1) 


The last two parameters to be discussed were introduced by General Electric and 
are called IDC max an d IDR^^. ** tlier e are five rings of total pressure measure- 
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ments in the plane of measurement, IDC max and IDR max can be defined by the follow 
ing expressions: 


IDC max = ^ 


jj- (IDCj + IDC 2 ), i (IDC 4 + IDC 5 )J 


IDR max = MAX ( IDR p IDR 5> 


where 


^ P t, av^ ~ ^ P t, min) T 

IDC-r = - 

P 

t, av, f 


P t, av, f “ ( p t, av) 

IDR = ± 

J P 

t, av, f 

and 

(P- ) average total pressure for ring J 

l,avj 

(P. min) minimum total pressure for ring J 

XX1JLI1 j 

P, f average of all total pressure measurements in the input measurement 

t, <IV, I 

plane 


The average total pressure for ring J, (P t ay )j,is stored in CMAVG(J); (P t min )j 
is stored in CMMIN(J); and P^ ay ^ is stored in TMAVG. The following FORTRAN IV 
program steps may be used to calculate IDC max and IDR max . For the sake of this 
example, we will assume that there are seven rings of input measurements and that the 
first and seventh rings consist of static pressure measurements. 


DO 10 J=2, 6 

RIDC(J)=(CMAVG(J) -CMMIN(J)) /TMAVG 
10 RIDR(J)=(TMAVG-CMAVG(J))/TMAVG 
RIDCXl=(RIDC(2)+RIDC(3))/2 .0 
RIDCX2=(RIDC(5)+RIDC(6))/2 . 0 
RIDCMX=AMAX1(RIDCX1, RIDCX2) 

RIDRMX=AMAX1(RIDR(2), RIDR(6)) 

A FORTRAN IV subroutine for these specific parameters may be seen in 
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appendix A The name ol the subprogram is DISPAR. The specific 

appendix a. MDnTTT ^hroue-h the labeled common block PARTRN . a 

transferred to the suprogr m ^ fQr user . wri tt en output statements to handle 

special section m LNPOUT available in the common 

the specific user-defined distortion parameters. Input data 

block DATTRN, but should not be altered. 
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appendix e 


TWO EXAMPLES OF INPUT AND OUTPUT 
Example 1 

This is an example of a circumferential total pressure distortion at the inlet of a 
jet engine. Figure 4 shows the instrumentation layout for the first example case. 
There is one ring of static -pressure measurements (at the outer duct wall), and there 
are five rings of total -pressure measurements. The pressure values corresponding to 
each measurement location are shown in table I. Ring 6 is actually a ring of static - 
pressure measurements. 

In this example, the interpolation process will continue through the center of the 
map (interpolation will not stop at the innermost ring), and (P. - P. „„ f )/P t „„ f (in 

1 l. 3.V . I l. 3.V , I 


Ring Radius, 
cm (in. ) 


1 8.00 (3.151 

2 16.03 (6.31) 

3 24.03 (9.461 

4 32.03(12.61) 

5 40.03(15.76) . 0° 



90° 


Figure 4. - Instrumentation layout for example 1. (Angles are measured clockwise 
from 0°. 
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Figure 6. - Output for example 1. 
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— EXAMPLE CASE NO. I — 

— WHICH SET TO 2, CENTER FILLED IN, AND DISTORTION PARAMETERS CALCULATED — 


CIRCUMFERENTIAL DISTORTION PARAMETERS 


ring no. 

I 

2 

3 

4 

5 

6 

7 

8 

9 

10 

thtmin 

195.69 

185.74 

172.56 

. 180.84 

182.23 

-0.00 

-0.00 

-0.00 

-0.00 

-0.00 

THTPLS 

164,31 

174.26 

187.44 

179.16 

177.77 

-0.00 

-0.00 

-0.00 

-0.00 

-0.00 

CMMIN 

45.300 

44.110 

42.410 

42.100 

42.070 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

CMMAX 

50.610 

52.390 

53.570 

51.730 

51.730 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

CMAVG 

47.563 

47.702 

47.874 

46.884 

46.867 

-0.000 

-0.000 

-0.000 

-0.000 

-0.000 

ICMMAX-CMMINI/CMMAX 

.10492 

.15805 

.20833 

.18616 

.18674 

.00000 

.00000 

.00000 

.00000 

.00000 

ICMMAX-CMMINI/CMAVG 

.11164 

.17358 

.23311 

.20540 

.20612 

.00000 

.00000 

.00000 

.00000 

.00000 

(CMAVG-CMMINI/CMAVG 

.04759 

.07529 

.11414 

.10204 

.10235 

.00000 

- . 00000 

.00000 

.00000 

.00000 





RADIAL DISTORTION PARAMETERS 


THETA POSITION ( OEG ) 

0 

45 

90 

135 

180 

225 

270 

315 

RMMIN 

42.025 

42.187 

42.270 

42.030 

41.625 

41.295 

41.280 

41.737 

RMMAX 

48.040 

46.472 

45.917 

45.362 

48.730 

52.607 

52.898 

51.017 

RMAVG 

46.365 

44.368 

43.300 

43.292 

46.507 

49.241 

49.786 

48.795 

(RMMAX-RMMIN) /RMMAX 

.12521 

.09219 

.07942 

.07345 

.14580 

.21504 

.21962 

.18189 

(RMMAX-RMMINI /RMAVG 

.12973 

.09656 

.08422 

.07696 

.15277 

.22974 

.23335 

.19017 

irmavg-rmmini/rmavg 

.59361 

.04914 

.02379 

.02914 

.10498 

.16136 

.17085 

.14463 


PARAMETERS PERTAINING TO THE 

ENTIRE ARRAY OF INPUT MEASUREMENTS 

TMM I N= 42.070 

TDP1=. 21467 

TMMAX- 53.570 

TDP2=. 24349 

TMAVG= 47.231 

T0P3=. 10927 


SMAVG= *1.808 

(c) Third page of output. 


— EXAMPLE CASE NO. 1 — 

— WHICH SET TO 2, CENTER FILLED IN, AND DISTORTION PARAMETERS CALCULATED — 


DEFINED DISTORTION PARAMETERS 


KD2 

1366.82698 


K-THETA 

5.78976 


I IOC > -MAX 

0.10143 


( IDRI-MAX 

5.00771 

(d) Fourth page of output. 



Figure 6. - Concluded. 
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percent) will be mapped. Simple and specific distortion parameters will be calculated. 
The specific distortion parameters are KDg, Kq, IDC max , and IDR m . 

The input for the first example is shown in figure 5, and the corresponding output 
is shown in figure 6. 


Example 2 

This is an example of a circumferential total temperature distortion at the inlet of 
a jet engine. Figure 7 shows the instrumentation layout for the second example case. 
In this example there are five rings of total temperature measurements. Since there 
is an inner boundary at this measuring station, the interpolation process will not con- 
tinue through the center of the map. In order to map the entire flow annulus, two 
imaginary rings (one at the outer wall and one at the inner wall) are generated. These 
fictitious measurement locations are shown as solid symbols in figure 7. 

The air temperature at the wall may be assumed to be at a value between the free- 
stream total and static temperatures . Therefore, each imaginary measurement will 
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TABLE II. - TEMPERATURE MEASUREMENTS FOR EXAMPLE 2 


Angle, 

Temperature, K (°R) 

deg 

Ring 1 

Ring 2 

Ring 3 

Ring 4 

Ring 5 

Ring 6 

Ring 7 

0 

328.1 (590.5) 

338.1 (608.5) 

337.7 (607.8) 

336.2 (605.2) 

335.1 (603.1) 

331.8 (597.3) 

321.8 (579.3) 

45 

322.7 (580.9) 

332.7 (598.9) 

331.9 (597.4) 

334.4 (601.9) 

334.4 (601.9) 

332.2 (597.9) 

322.2 (579.9) 

90 

295.8 (532.5) 

"30578 (550.5) 

299.4 (539.0). 

297 . 5 (535. 5) 

295.7 (532.3) 

295.4 (531.7) 

285.4 (513.7) 

135 

291.2 (524.2) 

301.2 (542.2) 

296.9 (534.5) 

295.9 (532.7) 

294.6 (530.2)' 

292. 7 (526. 9) 

282.7 (508.9) 

180 

291.4 (524.6) 

301.4 (542.6) 

297.3 (535.2) 

296.3 (533.3) 

297.6 (535.7) 

296.2 (533.2) 

286.2 (515.2) 

225 

292.3 (526.1) 

302.3 (544.1) 

298.4 (537.1) 

296.7 (534.0) 

295.1 (531.2) 

294.2 (529.6) 

284.2 (511.6) 

270 

291.0 (523.8) 

301.0 (541.8) 

298.0 (536.4) 

296.6 (533.8) 

295.4 (531.7) 

295.1 (531.1) 

285.1 (513.1) 

.315 

301.2 (542.2) 

311.2 (560.2) 

301.9 (543.4) 

298.3 (536.9) 

296.3 (533.3) 

294.5 (530.1) 

284.5 (512.1) 


FORTRAN STATEMENT 


■tt! 


ER r , NO PIS T O 


"SSSMI 


in 


Figure 8. - Input for example 2. 






















-- EXAMPLE CASE NO. 2 — 

— SPECIAL INPUT, OPEN CENTER, NO DISTORTION PARAMETERS -- 


RADIUS 

ANGULAR 

POSITION OF 

INPUT MEASUREMENTS, DEG 

16.00 

0 . 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 


21.10 

G • 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 


28.70 

0 . 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 


34.80 

0 . 

45. C 

90.0 

135.0 


225.0 

270.0 

315.0 


39.90 

0 . 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 


44.50 

0 . 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 


46. 50 

0 . 

45.0 

90.0 

135.0 


225.0 

270.0 

315.0 



RADIUS 

SPECIAL INPUT 




16.00 

328.10000 

322.70000 . 

295.80000 

291.20000 

291 


2 92. 30000 

291.00000 

301.20000 



21.10 

338 . ICoOO 

332.70000 

305.80000 

301.20003 

301, 


302.30000 

301.00000 

311.20000 



28.70 

337.70000 

331.90030 

299.49000 

296.90000 

297 


298.40000 

298.00000 

301.90000 



34. 80 

336.20000 

. 334.43000 

297.50000 

295.9000C 

296 


296.70000 

296.60000 

298.30000 



39.90 

335.10000 

334.40000 

295.70000 

294.60000 

297 


295.10000 

295.43000 

296.30000 



44. 50 

331.89030 

332.20000 

295.40000 

292.73000 

296 


294.20000 

295. 10300 

294.50000 



46.50 

321.80000 

322.20000 

285.43000 

282.70000 

286, 


284.20000 

285.10000 

284.50000 




(a) First page of output. 
Figure 9. - Output for example 2. 


180.0 

180.0 

180.0 

180.0 

180.0 

180.0 
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be given a value equal to 10 K (18° R) less than each corresponding adjacent total tem- 
perature measurement. For example, the imaginary measurement on ring 7 at 9 = 45 
(see fig. 7) would have a value equal to 10 K (18° R) less than the total temperature 
measurement on ring 6 at 0 = 45°. (Note that the program does not make this assump- 
tion. This is a user-defined assumption applied at the time of input data preparation. ) 
Rings 1 and 7 (the imaginary rings) should be treated as rings of static measurements 
when preparing the input. 

Table n shows the temperature values corresponding to each measurement loca- 
tion. A map of T +n(£T) will be generated, and £T is set to 3.0 K (5.4° R). 

The input for the second example is shown in figure 8, and the corresponding 
output is shown in figure 9. 


General Comments 

The lines of constant temperature in figure 10 were drawn by hand. The contour 
map used in figure 10 is the same contour map appearing in figure 6(a). Drawing in 
lines of constant mapping parameter by hand represents a useful way of looking at the 
contour maps. 

Questionable input measurements can be spotted by examining the contour maps. 
Referring to figure 6(b) (or fig. 10) near the 270° position, the small region of "5's" 
could be caused by a bad input measurement. Figure 6(a) shows that the region of "5's" 
is caused by the measurement of ring 4 (r=12.61 in.) at the 290° position. 

Note that the maximum "positive" pressure distortion (P. - P, f )/P, , in 

1 l, civ, I I, a-Vj I 

example 1 is over 13 percent and the maximum "negative" pressure distortion goes 
beyond -12 percent. For example 1, a value of RANGE = 20.0 was chosen. RANGE = 
20.0 was not large enough to cover the total spread of 25 percent occurring in example 1. 
Hence, in figure 6(b) areas where (P i - P^. av f )/P t av f is greater than 10 percent 
are represented by "9's" and areas where (P. - P, „ Tr ,)/P f f is less than -10 per- 
cent are represented by "J's". Also realize that making RANGE larger to cover the 
spread of percentage will reduce the resolution of the contour map. 
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Figure 10. - Example of hand-drawn contour lines. 
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APPENDIX F 


PROCEDURE WHICH ENABLES CALCULATION OF 
SIMPLE RADIAL DISTORTION PARAMETERS 

A geometric requirement of this program is that the input measurements have to 
be located on defineable rings. However, the measurements may be at any angular 
position on the rings. This means that a radial vector may exist at some angle 0 
without intersecting a single input measurement. (Remember that 0 is the reference 
angle for the input measurement plane and is measured clockwise from the 12 o’clock 
position.) 

Calculation of simple radial distortion parameters at given 0 positions requires 
interpolated values of the input measurement for each ring at each given 0 position. 
The circumferential interpolation procedure outlined in appendix B is used to find an 
interpolated value of the input measurement for each of the NR rings at eight discrete 
circumferential positions in the input measurement plane. These positions are 0°, 45°, 
90°, 135°, 180°, 225°, 270°, and 315° (Eight equally spaced 9 positions, 45° apart). 

Once interpolated values are known on each ring at each of the discrete 9 loca- 
tions, simple radial distortion parameters can be calculated for each discrete value of 
0. (The calculation of these simple radial distortion parameters involves only the 
values found by the circumferential interpolation process; no further radial interpolation 
is done.) The simple distortion parameters are shown in appendix D. 
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